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Force project monitors were Maj. K. B. Harwood (CF) and 
Mr. A. F. Money, both of the AEDC Deputy for Operations. 
The project sponsor was Mr. William Lane, AFFDL/FGC. 
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University of Tennessee Space Institute, Tullahoma, 
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1. INTRODUCTION 


Since the days of the early airships when the cal- 
culation of aerodynamic forces on the airships hulls was 
required, the methods for predicting slender body aero- 
dynamics have been extended in direct proportion to the 
increased performance of the various slender body air vehicle 
configurations. Increased maneuverability ard increased air 
speeds required the development of prediction techniques 
which were valid at higher angles of attack and supersonic 
Mach numbers. Continuing the trend to increased maneu- 
verability, applications for the flight of slender body 
configurations at angles to 180 degrees are being seen in the 
new generation flight vehicles currently being considered. 
The aerodynamic coefficient prediction techniques which were 
valid at angles of attack of approximately 20 or 30 degrees 
are no longer sufficient for preliminary design studies. The 
new concepts in slender body flight vehicles require new pre- 
diction techniques valid to 180 degrees angle of attack. 

An aerodynamic coefficient prediction technique 
which meets the above need is herein developed. The pre- 
diction technique relies on semi-empirical procedures for the 
calculation of normal force coefficients and pitching moment 
coefficients for slender body configurations with low aspect 


ratio fins. Additionally, the normal force coefficient for the 
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ee 


fin in the presence of the body is determined along with the 
longitudinal and lateral center of pressure of the force. 
Because the need to operate slender body flight 
vehicles at high angles of attack had not been established 
until recently, there was a complete dearth of experimental, 
parametric data for slender body configurations with low 
aspect ratio fins at high angles of attack. The small 
amount of data on configurations at high angles of attack that 
was available was highly configuration oriented and in some 
cases classified. The experimental, parametric data ob- 
tained as a part of this investigation is the first data of 
its kind and represents a broad base from which to develop 


empirical aerodynamic coefficient prediction methodology. The 


lack of parametric data for slender body configurations at high 


angles of attack before this time has hampered the development 
of empirical techniques for the prediction of aerodynamic 
coefficients. The future development of purely theoretical 
prediction techniques also requires a standard of comparison 
to evaluate the validity of their prediction capabilities. A 
data base which meets the requirements for developing empir- 
ical aerodynamic coefficient prediction techniques as well as 
providing a standard of comparison for theoretical techniques 


has been established by Baker qyt. 


lyumbers in parentheses refer to similarly numbered 
references in the bibliography. 
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In order to provide a zet of data which could readily be used 
to develop empirical methods for the calculation of the aero- 
dynamic coefficients, it was necessary to determine the effects 
of each major component of a generalized, finned, slender body 
on the forces and moments for the complete configuration. In 
order to determine these effects, a wind tunnel test was con- 
ducted on a series of parametric models. Provisions were made 
to test tail fins alone, using a reflection plane technique, 
so that the tail fin contribution to the total force and 
moment coefficient could be isolated. Also a generalized 
slender body was tested using a technique which allows the 
measurement of forces and moments on a single fin simultaneously 
with the total combined forces and moments acting on the com- 
plete configuration. The generalized configuration was 
tested both with and without tail fins. The change in the 
force and moment coefficients due to the addition of the fins, 
combined with the measured force and moment coefficients on 
the metric fin in the presence of the body and the fin alone, 
allows the determination of the mutual interference of the 
fins and body. 

A basic ogive cylinder body was designed along with 
a series of twelve tail fins which provided a parametric 
variation of the taper ratio, span ratio, and the aspect 
ratio of the fins. In order to provide data for a wide range 


of values of body length for the data set, tests were conducted, 
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using slender bodies with a total length of 10-, 12-, 12.66-, 
and i8-calibers (i.e., 10-, 12-, 12.66-, and 15-body diam- 
eters). Most of the data were obtained using the 10~caliber 
body and all 12 of the tail fins were tested with the 10- 
caliber body. The 12- and 12.66-caliber bodies were tested 
without tail fins and the 15-caliber body was tested alone 
and with the three tail fins having an aspect ratio of 0.5 
and with the fin having an aspect ratio of 2.0, span ratio of 


0.4, and taper ratio of 0.5. 
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II. EARLY HIGH ANGLE OF ATTACK PREDICTION METHODOLOGY 


The early requirements for the calculation of aero- 
dynamic forces and moments on slender body configurations were 
related to studies of rigid airships. A potential flow slender 
body theory was developed by Munk (2) which was quite adequate 
for angles of attack up to approximately five degrees. The 
potential flow equations for pitching moment and cross force 


per unit length, for airship hulls, are given in Equations 19 


and 23 of Reference (2). 


m = U2 —§—- (K, - K,) sin 2a (2.1) 
df = [v? - sin 2a a ax (2.2) 
Where 


p(K, - K,) is the apparent mass of the airship. 

Tsien (3) showed that Munk's formulation was valid 
for slender bodies at moderate supersonic speeds. 

At angles of attack above approximately five degrees, 
Munk's slender body theory begins to underpredict the normal 
force coefficient of a slender pody. Modifications to Munk's 
theory have been made to account for the viscous crossflow con- 
tributions to the forces and moments. This modification of 


Munk's theory, which has come to be known as the crossflow 
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theory, was developed in 1951 through the work of Allen and 
Perkins (4) wno assumed that the viscous contribution to the 
aerodynamic forces and moments on slender bodies et angie of 
attack could be separated from the potential flow contribu- 
tion. Thus, the force and moment equations were written as the 
sum of the potential term formulated by Munk (2) in 1924 and the 
viscous crossflow term formulated by Allen and Perkins (4). 
From the work of Ward (5), which showed that the potential 
cross force is directed midway between the normal to the axis 
of the body and the normal to the wind direction, Allen modi- 
fied the potential cross force term. Then by adding his 
viscous cross force term, Allen developed the following 


equations from Reference (4) for lift and pitching moment 


coefficients. 
Sb a £ 2 
Cc. == sin 20 cos xt Ca, + sin acos a (2.3) 


V- S) (2-x) 


Cm -[—e— | Sin 2a cos > 


s (X_ =X) 
p ™ sale 
+ Ca, $ =a Sin @ (2.4) 
Jorgensen (6! formulated the equations in terms of 


normal force, Cye and pitching moment, C_, for slender bodies 


m 
at high angles of attack and included the term n to modify 


the two dimensional crossflow drag, Ca, for the effects of a 
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? 
finite length body. The term n which modifies the two i 
dimensional drag coefficient to approximate the drag coeffi- 
cient for a finite length cylinder is determined from the data | 
obtained by Goldstein (7). Jorgensen's equation for normal i 
} 
force coefficient for the angle of attack range 0 £ a < 180 ‘ 
; degrees is given by 
- 
: s 3 'S ; i 
| The pitching moment equation for the angle of 
F ’ 
Q attack range 0 £ a £ 90 degrees is given by 
[ 
E Vv-S, (2 -X) , F 
k : poy en b m , a : : 
E i Cm sa ] stn (2a°) cos (5) 
4 
nS , 
b ; s (X,,-X) 
[ + ncq (=) —— sin” (a*) (2.6) 
F 
‘ i 
i and the pitching moment equation for the angle of attack ; 
range 90 < a <S 180 degrees is given by 
i f 
- . | V2 Sb%in} |, S a” : 
: Cr [==] sin (2a°) cos (3) j 
. ! 
: i 
és (2) (X.-%) in? (a’) (2.7) 
: +n sin a ° j 
—_ ™ : 
t 
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where 


a° =a 0 <a < 90 degrees 


a“ = 180 - a 90 <a < 180 degrees 
The location of the aerodynamic center given by 


c 
XCP = = (2.8) 5 
N 


Over the years the potential term of the equations 
has remained essentially unchanged and the greatest effort 


toward improving the cross flow theory has been associated 


with th- crossflow drag portion of the viscous term. The 
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most 1°:..2i. crossflow drag work is that of Fidler and 
Batemcn (8). The dramatic effect of the viscous contri- 
bution to the total forces and moments on a slender body is 
illustrated in Figure 1 taken from Reference (6). The 
variation with angle of attack of the first terms of 


Jorgensen's normal force equation (Equation 2.5) and pitching 
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moment equations (Equations 2.6 and 2.7) which represents 


Allen's modifications to Munk's slender body potential theory, 


7 
; 
7 
’ 
E: 
: 
’ 
; 
7 


are compared with the complete normal force and pitching 


moment equations. The shaded portion of the figure represents 


{ 


the viscous contribution to the crossflow theory. Because the 
potential term of the equation contains sin (2a) the term goes 
to zero at 90 degrees angle of attack, indicating that at 90 


degrees the coefficients are totally determined by the viscous 
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term. In the case cf the normal force, when the slender body 
is at 90 degrees angle of attack, the normal force is coinci- 


dent with the drag, thus the normal force is due to the skin 
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friction and separated flow over the lee side of the body, 
both viscous phenomena. 
Jorgensen's formulation of the crossflow theory 
equations for slender bodies at high angles of attack will be 
used as the basis for the body alone portion of the theory 
developed herein. 


Numerous calculation schemes nave been developed 
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for slender body configurations, both with and without fin or 
wings. Besides Jorgensen (6), others developing calculation 
procedures for slende~ boeies alone at. high angles of 

attack are Fidler and Batemen (8), Kellock and Miller (9), 


Thompson (10), and Gregoriou (11). Calculation procedures 
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Geveloped for slender bedies with fins and/cr wings 


i 
1 
{ 


have been written by Nikolitsch (12), Moore (13), Darling 

(14), baton (15), and Fidler and Batemen (16, 17), all of 

which are basically low to moderate angle of attack programs. 
High angle of attack programs for bodies with fins and/or wings 
were writcen by Saffell, et al. (18), Tipping, et al. (19) 

and Aiello (20). Mendenhall, et al. (21) provides a user's 
manual for four different computer programs for predicting 
aerodynamic coefficients. One program is the computer version 
of Jorgensen's method and the other three programs are low 

angle of attack body plus wing and fin programs. The details or 


three of the above high angle of attack aerodynamic coefficient 
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prediction programs will new be discussed. All three of the 
programs are based on various formulations of Allen's cross- 


flow theory. 


Lau? High Angie of Attack Aerodynamic Coefficient Prediction 
Frojran 


The LSU high angle of attack aerodynamic coefficient 
prediction program (Tabie 1A), was written by Captain Robert 
E. Kellock (9). The program calculates aerodynamic coeffici- 
ents for a body alone configuration. The primary calculation 
scheme was formulated by Kelly @2). The normal force is 
assumed to vary along the body and the total normal torce and 
pitching moment are determined by integrating the normal force 
distribution along the body. The body is divided into elements 
and the crossflow drag is assumed to vary from one element to 
another. The crossflow drag of each increment is calculated 
as the drag of an impulsively started cylinder. The time 
parameter associated with the impulsively started cylinder is 
related to the axial velocity and the axial position of the 
incremental element. The drag data of Schwabe (23) for an 
impulsively started cylinder were used to modify the steady- 
state crossflow drag coefficient, taken from Hoerner (24), as 
a function of crossflow Mach number and crossflow Reynolds 


number. The crossflow drag at each element is used with the 


érsu is an abbreviation for Louisiana State 


University. 
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crossfiow theory to calculate a distribution of normal force 
along the body. The three-dimensional effects of the finite 
length cylinder are accounted for by the values of n given by 
Goldstein (7). The predictions of normal force and pitching 
moment made with the LSU program were compared with the data 


from the high angle of attack data set in Reference (1). ‘ 
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NSRDC” High Angle of Attack Aerodynamic Coefficient j 
Prediction Program ae | 


The NSRDC high angle of attack aerodynamic coefficient 


ogTrecime erie 


prediction program (Table 1B) was written by Saffell, Howard, ' 
and Brooks (18). It consists of a main program and three ' 
subroutines. The main program calculates the body alone lift, ! 
drag and pitching moment coefficients using the crossflow drag 
theory. Three-dimensional effects are accounted for dy the 
values of n given by Goldstein. The values of Ca, are deter- 
mined from Allen and Perkins (25). The lift coefficients for 
the wings and tails and the wing-body and tail-body inter- 
ference factors are calculated using linear relations for the 
lifting surfiaces modified by a function of sin a and con- 
strained tc be zero at ov = 90 degrees. The lift of the tail 
surfaces is further modified to account for the vortices shed 
from the wing surfaces. The pitching moment calculations are 


based-on the formulation of the equation by Allen. The drag 


3ysRDC is an abbreviation for Naval Ship Research 
and Development Center. 
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calculations are based on the methods from the USAF 

DATCOM (26). Calculations were made using the NSRDC Program 
and compared with the data from the high alpha data set in 
Reference (1). The program computes the static aerodynamic 
coefficients for bodies with wings, tails, and strakes, 
including control surface deflections, for an angle-of-attack 


range fror -180 to 180 degrees. 


Computer Aided Missile Synthesis Program (CAMS) 


The CAMS program written by Tipping, et al. (19) is 
used to design a complete missile system. The program uses an 
iteration scheme to vary the many design parameters associated 
with a missile design. The aerodynamics module of the CAMS 
program (Table 1C) can be used by itself as a coefficient pre- 
diction program. It is the most complex of the coefficient 
prediction programs discussed here because of the number of 
different configurations that can be analyzed. The program is 
similar to the other programs in that it requires geometric 
characteristics of the model and flight conditions as program 
inputs to calculate the aerodynamic coefficients. The cal- 
culation scheme is similar to the other prediction programs 
in that the potential and viscous crossflow terms are com- 
bined linearly for the body alone normal force and pitching 
moment coefficients. Empirical data are used extensively 


throughout the program for determining the effects of fins, 
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strakes, etc. For each Mach number range, transonic, super- ; 
sonic, and hypersonic, different relations are used for the j 
“ 

linear and nonlinear iift contributions of the various com- : 


ponents of the configuration. For the body aione, the cross- 
flow theory is used with modification for boattail effects. 


Calculations were made for the high angle of attack data 


base models and compsred with the experimentally determined 


coefficients in Reference (1). 
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IIIT. EXPERIMENTAL TECHNIQUES 


Wind Tunnel vescription 

The tests whicr established the data base were 
conducted in the Aerodynamic Wind Tunnel (4T) of the 
Propulsion Wind Tunnel Facility (PWT) and the Supersonic 
Wind Tunnel (A) of the Von Karman Gas Dynamics Facility 
(VKF)} at the Arnold Engineering Development Center (AEDC). 

The Aerodynamic Wind Tunnel (4T) is a continuous 
flow, closed-circuit, variable density wind tunnel capable 
of operating at stagnation pressures of 300 to 3700 psfa 
and at Mach numbers from 0.1 to 2.0. The Mach number is 
continuously variable from M = 0.1 to 1.3 and nozzle 
inserts can be installed to reach M = 1.6 and 2.0. The 
nozzle consists of a contraction section which serves as 
a transition from a circular stilling chamber to a 
rectangular nozzle. The solid block, sonic nozzle is 
composed of flat sidewalls and contoured top and bottom 
walls. The desired Mach number is generated by controlling 
the pressure ratio across the nozzle and by regulating the 
plenum suction. The test section is 4 ft square and 12.5 ft 
long. It is equipped with four variable porosity walls 
adjustable from 0 to 10%. The two test section sidewalls 
are fixed and the top and bottom walls are adjustable +1/2 
degree from parallel. The test section is enclosed by a 14 x 14 


ft square plenum chamber which can be evacuated allowing part 
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of the tunnel main flow to be removed through the test 
section walls to both generate supersonic flow and reduce 


wall interference. 


The Supersonic Wind Tunnel (A) is a continuous, 
closed-circuit, variable-density wind tunnel with an 
automatically driven flexible-plate-type nozzle and a 
40 x 40 in. test section. The tunnel can be operated at 
Mach numbers from 1.5 to 6.0 at maximui stagnation pres- 
sures from 4,200 to 28,000 psfa, respectively, and 
stagnation temperatures up to 750°R (M = 6). Minimum 
operation pressures range from about one-tenth to one- 
twentieth of the maximum at each Mach number. Although 
Tunnel A is primarily a supersonic tunnel it can be 
operated subsonically from Mach numbers 0.2 to 0.8 by 
opening the throat (M = 1.2 setting) and closing the 
diffuser until the tunnel chokes at that point. The 
tunnel is equipped with a model injection system which 
allows removal of the model from the test section for 


model changes while the tunnel remains in operation. 
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Wind Tunnel Models 
Two different type models were tested in Tunnels 4T and A. 
The first type of model was a strut supported siender body 
tested both with and without tail fins. Body alone, body 
plus fin and fin the presence of the body data were obtained 
with this model. This model was designed so that it had less 
than 1% blockage in the tunnel at 90 degrees angle of attack. 
The second model type was a reflection plane mounted fin 
which was tested to obtain fin alone data. The fins had less 
than 0.6% blockage at 90 degrees angle of attack and the total 
blockage including the reflection plane was less than 1.58%. 
Slender bodies having total lengths of 10~, 12-, 
12.66-, and 15-calibers were tested. Each body consists of a 
2.5-caliber sharp tangent ogive nose, designated N2, and a 
cylindrical afterbody. The basic dimensiones 0° the ogive 
nose and cylindrical afterbodies are given in Figure 2. The 
7.5-, 9.5-, 10.16-, and 12.5-caliber afterbodies are desig- 
nated as Al, A2, A4, and AS, respectively. 
Two different size tail fins were used in the test. 
The tail fins tested with the reflection plane were geomet- 
rically similar to but three times larger than the fins 
tested with the slendex body. The dimensions of the fins 


used with the reflection plane are given in Table 2. The 
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4 
dimensions of the fins tested in conjunction with the slender 
body are shown in Figure 3. Because of the smaller size, the | 
fins tested with the slender body are not as detailed as 


those tested on the reflection plane. The tail fin configura- 


en oe ed 


tions are identified by the designation Txx, where the two 4 
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chord and were located at 55 and 62% of the root chord, 


digit number (xx) is given in Table 2 and Figure 2? along with { 

the fin dimensions. Thus, the total configuration designation i 

for a typical configuration is N2A1T32. The nomenclature de- 

scribed above was adopted to be consistent with that for a ; 

similar set of models tested extensively at relatively lower 

angles of attack by Fidler and Bateman (8): § 

i For each fin type, rectangular, trapezoidal, or E | 
triangular, the hingeline, HL, was located at a different : 
; position. For the rectangular fins, represented by A = 1.0, 
the hingeline was normal to the root chord and located at ; 
45% of the root chord measured from the leading edge. For 7 
i the trapezoidal fins and the triangular fins, \ = 0.5 and 
: 4 = 0.0, respectively, the hingelines were normal to root i 
| 


respectively. 


Reflection Plane Model Installation 
The installation of the reflection plane on the 
: main sting support system in the test section of Tunnel 4T ‘ 


and Tunnel A is shown in Figure 4. The details of the 
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reflection plane assembly ére shown in Figure 5. To permit 


whe 


testing through the desired angle-of-attack range, the 
reflection plane assembly includes a drive mechanism which 
allows the initial fin angle relative to the reflection plane 


to be indexed remotely. The assembly also contains provisions 
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for indicating the discrete initial fin angles of 0, 15, 30, 
60, 90, 120, 150, and 180 degrees. After remotely setting 
the initial fin angle, the angle-of-attack sweep is made, 


using the main model support system. 


Generalized Slender Body Installation 


Two strut mounting techniques were used to support 
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the generalized slender body models and achieve the angle-of- 
attack range from 0 to 180 degrees. A typical installation 
cf each mounting technique is shown in Figure 6 for the 

&£/d = 10 configuration in Tunnels A and 4T and Figure 7 for 


the £/d = 15 configuration in Tunnel 4T. The details of 
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the model installations are shown in Figures 8, 9, 10, and 11 


for %/d = 10, 12, 12.66, and 15 models, respectively. The 


Poe 


model support system includes a clutch face joint which 
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allows the initial or prebend angle of the model to be varied 
in five-degree increments from 2.5 to 177.5 degrees. 


Instrumentation 


pe eck apie oom 


Aerodynamic loads on the fins during the reflection 
plane portion of the tests were measured, using a three- 
component internal strain-gage balance. The positive direc- 


tion of the forces and moments are shown in Figure 12. 
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Strain gages were attached to the sting for the transonic 
tests so that the angle of attack of the fins could be cor- 
rected for both balance deflections and deflection of the 
sting caused by loads on the non-metric portion of the re- 
flection plane and sting. 

Aerodynamic loads on the slender body plus fin 


models were measured with a six-component internal strain- 


gage balance. In addition, a three-component internal strain- 


gage balance, mounted at the rear of the model, measured the 
aerodynamic loads on one fin. The positive directions of the 
forces and moments are shown in Figure 13 for the slender 


body model. 


Test Procedure 

The reflection plane tests were conducted in two 
phases: the first at transonic speeds in Tunnel 4T and the 
second at supersonic speeds in Tunnel A. In both tunnels the 
initial angle of the fin was indexed remotely in 30-deg in- 
crements, and at each setting an angle-of-attack sweep was 
made using the automatic pitch support system of the tunnel. 
The pitch-pause technique, in which the automatic pitch 
mechanism stops at each discrete angle of attack to record 
data, was used in Tunnel 4T. The continuous sweep technique, 
in which the automatic pitch mechanism moves continuously 
while the data are recorded, was used in Tunnel A. Several 


pitch-pause sweeps were made to check the continuous sweep 
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data. In Tunnel 4T the data were recorded while the angle 
of attack was increasing, whereas in Tunnel A the data were 
recorded while the angle of attack was decreasing. The bound- 
ary layer on the fins was allowed to transition naturaily. 

The generalized missile models were tested in a 
manner very similar to the reflection plane-mounted fin 
models. The transonic tests were conducted in Tunnel 4T, and 
che supersonic tests were conducted in Tunnel A. The pitch- 
pause technique, with the model pitching in the positive di- 
rection, was used in Tunnel 4T. The continuous sweep tech- 
nique, with the model pitching in the negative direction, was 
used in Tunnel A. The prebend angle (a,5) was manually set. 
by adjusting the clutch face joint to the desired angle. 
During the tests in Tunnel 4T, the prebend angles were ad- 
justed in 30-deg increments. The angle of attack was varied 
during the test using the main support pitch system. The 
35-deg movement of the Tunnel 4T pitch sector allowed a 5-deg 
Overlap in the data at each prebend setting. During the 
tests in Tunnel A, the prebend angles were adjusted in 25-deg 
increments. The angle of attack was varied during the test 
using the main support system. The 20-deg movement of the 
Tunnel A pitch sector left a 5-deg gap in the data at each 
prebend setting. 

Data were taken with artificially induced transi- 


tion on the body. The boundary-layer trips consisted of two 


32 


a AR RANA Re tt oe DAE GUE BE eT CR 


Tyee iste aga o> 


a NI aa PIU RTE St PCY EO SM 


w 


longitudinal rays of No. 180 grit applied in 0.1 inch wide 
strips itocated 30 deg either side of the windward ray (Fig- 
ure 14). The grit was applied in longitudinal rays because 
at high angles of attack the conventional application of 
grit, a small ring of grit around the nose near the tip, 
would be ineffective. The size of the grit to be used was 
determined by the method of Braslow and Knox (27). During 
the test, other grit patterns and sizes were investigated 
with no appreciable difference in the data from that obtained 
with the two longitudinal rays of grit. The presence of the 
grit caused a noticeable difference in the data at Mach num- 
bers 0.6 and 0.8 but no noticeable difference at higher Mach 
numbers. 

The grit which was applied to the model to artifi- 
cially induce transition had to be reapplied at Frequent in- 
tervals. The conventional technique of applying grit with 
Polaroid® print coater as the adhesive proved inadequate, 
causing the grit to be blown from the model by the air stream. 
Numercus adhesives were tried with Eastman 910° adhesive be- 
ing the most effective. The area where grit was applied was 
outlined with masking tape and the adhesive was applied to 
the surface. The grit was then blown onto the surface and the 
masking tape was removed. A uniform distribution of grit 


with the individual particles not touching each other was 


strived for. 
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Test Conditions 

The fin alone tests, using the reflection plane 
technique, were conducted at Reynolds numbers which varied 
with Mach number. The unit Reynolds number varied from 
1.1 x 10 at M = 0.6 to 3.2 x 10° at M = 3.0. It was planned 
to conduct the fin alone tests such that the Reynolds number 
based on fin chord at each Mach number was approximately the 
same as the Reynolds number based on fin chord in the body 
plus fin tests. This match was possible at Mach numbers 0.6, 
0.8, and 0.9; however, due to operating restrictions, it was 
necessary to conduct the tests at higher Mach numbers at 
higher Reynolds numbers. 

The tests on generalized slender body configurations 
both with and without tail fins were conducted at Mach num- 
bers of 0.6, 0.8, 0.9, 1.0, 1.15, 1.3, 1.5, 2.0, 2.5, and 3.0. 
For most test conditions, a unit Reynolds number of 4 x 10° 
was maintained resulting in a Reynolds number based on body 
diameter of 4.2 x 10°. For some of the tests on configura- 
tions with large fins, the unit Reynolds number was reduced 
to 1.5 x 10° to prevent overloading the fin balance, result- 
ing in a Reynoids number based on body diameter of 2.6 x 10°, 
For selected configurations, the unit Reynolds number was 
varied from 2 x 10 to 4 x 106 providing a variation of 
Reynolds number based on body diameter from 2.1 x 10° to 
4.2 x 10°. 
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The model was inspected during every model change 
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and the longitudinal rays of transition grit on the slender 
body model had to be reapplied approximately every two hours 
of tunnel operation. 

The total temperature in the tunnel was maintained 


at approximately 100 to 110 degrees F, except for a few cases 
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where the temperature had to be raised as high as 130 degrees 
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F to prevent the formation of fog at supersonic Mach numbers. 


The teste were conducted with no visible moisture in the test 


section. 
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Precision of Measurements 
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Uncertainties (bands which include 95% of the cali- } 
bration data) of the basic tunnel parameters (P,) and (M) 
were estimated from repeat calibrations of the instrumenta- . 
tion and from che repeatability and uniformity of the test 
section flow during tunnelcalibrations. The uncertainties 
were combined using the Taylor series method of error propaga- 
tion to determine the precision of the reduced parameters 


presented in Tables 3 and 4. 


35 


? 
& 
é 
i 
i 
i 
t 
& 
i 
E 
i 
| 


Table 3 
Fin-Alone Data Precision 
0.1619 |; 0.128 
0.1090 | 0.086 
0.0960 | 0.075 
0.0600 | 0.048 
0.0480 | 0.036 
0.0430 | 0.035 
0.0018 | 0.0013 
0.0019 | 0.0014 | 
| 36 
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0.046 
0.032 
0.028 
0.017 
0.013 
0.013 
0.0095 


0.0095 
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IV. AEDC HIGH ANGLE OF ATTACK DATA BASE 


Most of the high angle of attack data which 
were available prior to this investigation were primarily 
for supersonic Mach numbers and for body alone type con- 
figurations as in References (28) through (40}. The few 
examples of data for finned slender bodies are given in 
References (41) through (47). Again, most of the data 
are for supersonic Mach numbers. References (44), (45), 
and (46) contain data at transonic Mach numbers but all 
three tests were highly configuration oriented and had 
either large aspect ratio tail fins or ringtails. Only 
keference (47) contains transonic data for a slender body 
with low aspect ratio tail fins. 

Because of the complete lack of parametric data at 
high angles of attack, it was necessary to conduct an exten- 
sive wind tunnel test program to establish a set of data from 
which to develop an empirical aerodynamic coefficient pre- 
diction technique. It was desired to have a wide range of 
parameters typical of finned slender bodies. Thus the main 
parameters varied in the wind tunnel program were selected to 
be the total body length, the ratio of the body diameter ty 
the total fin span, the taper ratio of the fins and the aspect 


ratio of the fins. A basic ogive-cylinder configuration was 
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selected as the slender body to be tested. Mounted on the 
slender body was a set of four fins in a cruciform plus 
orientation. The basic configurations were tested at Mach 
numbers over the range from 0.6 to 3.0 and at angles of attack 
from 0 to 180 degrees. The three different types of data 
that were obtained in the data set, fin alone, body alone, and 
body plus fin, will be discussed individually. 
Fin Alone Data 

The fin alone data measured for the high angle of 
attack data set consisted of the fin normal force, hinge 
moment and root bending moment. Dividing the moments by the 
normal force resulted in the longitudinal and lateral centers 
of pressure of the normal force. Typical data obtained for 


three fin alone configurations are shown in Figures 15, i6, 


oe 


and 17 at Mach numbers 0.8, 1.3 and 2.0, respectively. The 
fins have an aspect ratio of 2.0 and semi-span as indicated 

in Table 2. Subsonically the sudden decrease in fin normal 
force at approximately 20 to 30 degrees angle of attack 
represents the typical stall condition with the triangular, 

A = 0.0, fin reaching the highest fin normal force before tie 
stali and the square, 4 = 1.0, fin achieving the lowest. 

Tests with the aspect ratio 0.5 fin showed the opposite trend. 
This trend with aspect ratio is reasonable since Bradley, et al. 
(48) using Polhamus' leading edge suction analogy (49) showed 


that while the potential lift of a rectangular fin decreased 
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\ with decreasing aspect ratio the tip vortex lift increased 
greatly providing an overall increase in lift with decreasing 
aspect ratio. Polhamuse (49), however, showed that the vortex 


lift of a triangular fin remained relatively constant with 


decreasing aspect ratio while the potential lift decreased 
with decreasing aspect ratio resulting in an overall decrease 
in lift for triangular fins with decreasing aspect ratio. 
After the stall, the normal force of all three fins increases 
to a maximum at 90 degrees angle of attack. As can be seen, 


the stall occurred at higher angles of attack for decreasing 
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taper ratio, a trend also predictable by the leading edge 
} suction analogy. : ae 
: As the angle of attack of the fins is increased 
beyond 90 degrees, the subsonic weal again appears, but. since 
all three fins have straight trailing edges, they all behave 
similar to a A = 1.0 fin with the stall occurring at about 
180 degrees minus the square fin stall gngve: For the tests 
conducted with fins of smaller aspect ratio, it was noted that 
the subsonic stail occurred at increasing angle of attack for 
decreasing aspect ratio. This trend is discussed by Bradley, 
et al. 

At low supersonic Mach numbers, the typical stall does 
not occur and the normal force increases smoothly to the maximum ‘ 
at $0 degrees angle of attack. A gradual increase in maximum 


fin normal force was nocred with increasing Mach number, up 


to M= 1.3. The maximum fin normel force was approximately the 


same level at M = 1.15 and 1.3. 


40 


3 
ra 
& 
z 
¢ 
3 
: 
¥ 
i 
2 


t 
ie i al tee ca, pathic ak 


5 ee OR TA a TN TE a te UE a Siac. VEE oe ewe) x Sips FES Bhs aot * 


~ 


ai 


{ 

| At higher supersonic Mach numbers, the nermal 

7 

i force as a function of angle of attack began to decrease at 
angles of attack between 50 and 130 degrees with the 


: minimum occurring at 90 degrees. This dip in the normal 


selec bat ed Ete 


Ft force can be seen in Figure 17 and is the result of 
separated flow on the reflection plane ahead of the fin. 


As the angle of attack of the fin reached approximately 50 
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degrees the flow on the reflection plane began to separate. 
The separated flow increased with increasing angle of attack 
of the fin, Figure 18, until the maximum region of separation 
was reached with the fin at% = 90 degrees. The separated 
region decreased in size as the angle of attack was increased 
; beyond 90 degrees until the separation disappeared at 
approximately 130 degrees angle of attack. The measured 
normal force at 90 degrees angle of attack decreased with 
increasing Mach number. The data affected by the separation 
were corrected before they wee used in the development. of the 
prediction technique. The details of this correction are 
: given in the next section. 
; The variation of root bending moment with angle of 
attack for all of the tin alone configurations followed very 
t closely the fin normal force variation with the maximum value 
: decreasing with decreasing taper ratio. The resulting lateral 


center of nressure occurred at approximately the lateral cen- 
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troid of the fin. The longitudinal center of pressure on the 


fins began at a forward position on the fin and moved aft with 
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i increasing angle of attack to the approximate longitudinal cen- 


troid of the fin at 90 degrees angle of attack and then con- 
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tinued on toward the trailing edge as the angle orf attack was 


increased beyond 90 degrees. The trend in longitudinal and 


oe nt ae a cate, 


lateral centers of pressvre was the same for all aspect 


ratios and Mach numbers. 
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Corrections to Fin Alone Data 
As mentioned earlier, separated flow on the reflection 
plane resulted in erroneous data for the fin normal force and 


the root bending moment at supersonic Mach numbers (1.5 <M < 3.0). 


From oil flow movies it was determined that significant separation 
was present. between the angles of attack of 50 and 130 degrees. 
Therefore, the data at angles of attack lese than 50 degrees 

and more then 130 degrees were assumed to be correct. ‘he 
correction to the fin normal force consists of determining a 

value of fin normal force for each fin at an angle of attack 


of 90 degrees and fairing a smooth curve from the correct 
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data at 50 degrees through the determined value at 90 degrees 


to the correct data at 130 degrees. The y-center of pressure 
of the fin normal force was corrected in the same manner. 
No effect was noted in the x~center of pressure of the fin 


i normal force due to separation. Thereforeno correction was 


made. 
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} An example of the corrected data is shown in Figure 
fs 19 where the fin alone data before and after the modification 


are compared with the installed fin data. 
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The maximum value of fin normal force, which was | ! 
: assumed to occur at 90 degrees angle of attack was determined 
from the calculation of a flat plate at 90 degrees, reference 
(24) and from a set of data for delta wings alone, Falunin, 
' 4 et al. Reference (50). The calculated value for a flat 
| plate was CNe = 1.7 over the Mach range from 1.5 to 3.0. | 
t The delta wing data, for an AR = 0.706 wing had values of Cy ; i 
of 1.38, 1.754, 1.759 and 1.694 at Mach numbers of 1.5, 2.0, i 
i 2.5, and 3.0, respectively. Figure 20 shows data from the high 
§ angle of attack data set for a series of AR = 1.0 fins { | 
. compared with the flat plate calculation and the delta wing data. ‘ 
Since the data from the high angle of attack data set indicate 
| that there is a decrease in Cy as a function of taper ratio, : 
; the assumed values of fin coeead force indicated by the 
: solid symbols are 1.7 for the rectangular fin, 1.65 for the 
i trapazodial fin and 1.6 for the triangular fin. These 
f values are assumed constant with Mach number from Mach numbers 
1.5 to 3.0. The assumed values of cNe compared with other data : 

i from the high angle of attack data set are shown in Figures 


20b and c for aspect ratio 2.0 and 0.5, respectively. 
The correction to the y-center of pressure was 
made over the same angle of attack range that the normal force 


correction was made. In order to determine the y center of 
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pressure at 90 degrees angle of attack, the assumption 


ON PIEPER EN SND Mtn cag SS Pen Ten 


was made that the center of pressure was 10% inboard of the 


centroid of the exposed fin area. A smooth curve was then 
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assumed beginning with the measured data at 50 degrees 
angle of attack, passing through the assumed point at 

90 degrees and continuing to the measured data at 

130 degrees angle of attack. An example of the 

corrected data compared with the fin alone and installed fin 


data is shown in Figure 21. 
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As has been stated, no correction was made for 


the x center of pressure. A comparison of the x center of 


pressure measured for the fin alone and installed fin cases 


is shown in Figure 22. 


Body Alone Data 

The body alone data measured for the high angle 
of attack data base consisted of normal force, pitching 
moment, side force, yawing moment, rolling moment and 
axial force. Only the pressure of the normal force was 


determined by dividing the measured pitching moment by the 
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measured normal force. Typical examples of the body alone 

data are shown for four different length models in Figures 23 
and 24 for Mach numbers 0.8 and 1.3, respectively. The magnitude 
of the normal force at both subsonic and low supersonic Mach 
numbers increases with body length approximately in proportion 

to the increase in planform area associated with the increase 

in body length. At subsonic Mach numbers, the normal force 


increases smoothly up to approximately 60 degrees angle of 
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attack. In the angle of attack range from 0 to 40 degrees, 
Thompson (10) has described the wake behind a slender body to 
be steady while the wake in the angle of attack range from 40 
to 60 degrees is describec as quasi-steady. In beth the 


Steady and quasi-steady regions the normal force are 
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expected to increase smoothly with angle of attack. 
Above 60 degrees, Thompson describes the wake as unsteady 


and indicates that the level of normal force should decrease. 
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The data in Figure 23 for the different length bodies do not 


decrease above 60 degrees but it does level off with a 
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gradval increase to the maximum value at 90 degrees. 
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Also noted in Figure 23 for the 2£/d = 10 configuration, 
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N2A1T00, the data obtained at 90 degrees angle of attack 
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with the aft mounted strut, Figure 8, does not agree with 


Po ae 


the data obtained using the nose mounted strut, Figure 
9. This mismatch in data which occurred at Mach numbers 0.6 


and 0.8 is an indication of support interference at subsonic 
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Mec’ numbers. The support interference problem will be 


discussed in a later section. 


i 
| 


At supersonic Mach numbers, the normal force in- 


ATES 


; creases smoothly to the maximum at 90 degrees angle of attack 
: for all of the body lengths tested. The maximum normal force 
for each ry length increased gradually with increasing sub- 
sonic Man number. The maximum normal force was apprceximately 
; the same for Mach numbers 1.0 and 1.15 and began a gradual 


decrease it vej. with increasing supersonic Mach numbers. 
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No indication of support interference was observed at super- 
sonic Mach numbers. 

The body alone pitching moment showed the same 
trend, with body length as the normal force, with the maximum 
and minimum levels increasing with body length in approximate 
proportion to the increase in body planform area fore and aft 
of the moment reference center. The maximum and minimum 
levels of pitching moment for each body length increased in 
magnitude for increasing Mach number up to M= 0.9. At Mach 
number 1.0, the magnitude of the maximum value of pitching 
moment began to decrease with increasing Mach number. Beyond 
M = 1.0, the magnitude of both the maximum and minimum values 
of pitching moment showed a decrease with increasing Mach 
number. The angle of attack at which the maximum pitching 
moment occurs decreases slightly with increasing Mach number 
while the angle of attack at which the minimum value of 
pitching moment occurs increases slightly with increasing 
Mach number. 

Interference in the form of mismatch in the data 
obtained using the forward and aft struts was observed in the 
pitching moment at M = 0.6 and angle of attack of 90 degrees. 
This interference is attributed to the proximity of the model 
base to the wind tunnel wall for the forward strut supported 
£/d = 12.66 and 15 slender body models, Figures 10 and 11. 


The position of the £/d = 15 model relative to the wind 
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tunnel wall at 90 degrees angle of attack is seen in Figure 


7. At a Mach number of 0.8 the mismatch was almost undec- 


table, Figure 23, and at higher Mach numbers the data from 
the forward and aft strut mounted models were in excellent 
agreement at 90 degrees. 

Throughout the data base, excellent matching was 
obtained at the angles of attack where a strut or prebend 
angle change was made and data were obtained at overlapping 
angles of attack. In only a very few cases was mismatch 


observed. For every slender body configuration tested over 
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the complete angle of attack range from 0 to 180 degrees, 


four prebend angle changes were made and one strut change was 
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made. Thus five overlaps cccur for each configuration tested 
over the complete angle of attack range. The excellent match- 


ing was obtained in all coefficients measured. 


Body Plus Fin and Installed Fin Data 
The body plus fin data measured for the high angle 
of attack data base are the same as described for the body 
alone case. Again, only the normal force and pitching moment 
data will be discussed here. Body plus fin data typical of 
the data in the data base are shown in Figures 25, 26, and 27 
Mach numbers 0.8, 1.3, and 2.0. The data shown are for three con- 
é figurations, each having an £/d = 10 body and tail fins hav- 
ing taper ratios of 0.0, 0.5, and 1.0. Each tail fin has an 


aspect ratio of 2.0 and a span ratio of 0.4. As can be seen 


47 


Dt Cec a a ee 


Lae ners non tt arememsae bent tibe £28 i seed! Pipa clabl Kt dette Beso Wetieieli, cane ARE ty Se 2: nf o> 6.22 . ~ 


ace ai! eater sain eine ode te al cai ed a iE i tcl De hl te tlle La ill re bi I ea Oe 


j 
j 
| 
} 
{ 
{ 
, 


ern a eestenammaeat 


| 


$0. Vee REPO Centnstematmereememr cere enn a ne te ne mee ee 4 
( 
i 
4 


from the figures, the body plus fin normal force data are 
essentially independent of taper ratio. This trend is typical 
of most of the data except for the data obtained with the 

AR = 2.0 and d/b' = 0.3 fins, where at Mach numbers 0.6 and 
0.8 significant vortex lift was developed on the triangular 
fin at angles of attack up to the fin stall angle causing a 
significantly higher fin normal force resulting in a higher 
total normal force. It should be noted that the semi-span of 


the d/b' = 0.3 fin is greater than for the other fins. Thus 


the fin protruded further into the freestream, resulting in the 
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least body influence. For the triangular fins with smaller 
semi-spans, the presence of the body appears to decrease the 
amount of vortex lift developed. 

For fins having a constant aspect ratio of 2.0 and 


a constant taper ratio of either 0.0, 0.5, or 1.0, the maximum 
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total normal force was increased by decreasing the span ratio 
from 0.4 to 0.3. Decreasing the span ratio physically means 
increasing the semi-span; therefore, if the aspect ratio is 


held constant, the smaller span ratio results in a larger area 
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: fin. The increased force with decreasing span ratio then is 
probably caused both by increasing the fin area and by moving 
the centroid of the fin further out into the airstream. 

The data obtained for the configurations having a 


constant span ratio of 0.5 and constant taper ratios of either 
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0.0, 0.5, or J)..0 show that decreasing the aspect ratio from 
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1.0 to 0.5 resulted in an increase in maximum total normal 
force. Just as above, decreasing the aspect ratio with con- 
stant semi-span results in an naceeae® in fin area and a 
resulting increase in total normal force. 

For configurations having the fin area aporoxi- 
mately constant, changing the span ratio, taper ratio, and 
aspect ratio had relatively little effect. The effect of Mach 
number and angle of attack on the body plus fin maximum total 
normal force is essentially the same as for the body alone 
configuration. 

The pitching moment unlike the total normal force 
has a slight dependence on taper ratio at subsonic Mach nun- 
bers for all of the taper ratio 0.0, triangular fins. For 
the aspect ratio 2.0 fin with span ratio 0.3, the increase in 
vortex lift on the triangular fin resulted in a significantly 
more negative pitching moment up to slightly above the fin 
stall angle. The aspect ratio 2.0, span ratio 0.4 fin, and 
the aspect ratio 1.0 fins have pitching moments siightly more 
negative for the triangular fin due to the longitudinal cen- 
troid of the fin being further aft than for the other two 
fins. The aspect ratio 0.5 triangular fin, just as in the 
fin alone case, had a lower value of fin normal force than 
the taper ratio 1.0 and 0.5 fins. Therefore, the resulting 
pitching moment had a less negative maguitude than the other 


two fins. As would be expected, changing the fin area has a 
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significant effect on the pitching moment with largest fins 
having the most negative pitching moment. 

The pitching moment as a function of angle of 
attack, as seen in Figure 25, smoothly decreases to a minimum 
at approximately 35 to 40 degrees angle of attack. The 
pitching moment then increases to a maximum at approximately 
55 degrees angie of attack followed by another decrease to a 
second minimum of approximately 120 degrees. This reversai 
of slope of the pitching moment at 35 to 40 degrees angle of 
attack diminishes with increasing Mach number. At a Mach 
number of 1.15 an inflection point in the curve occurs at 35 
to 40 degrees angle of attack and at supersonic Mach numbers 
the pitching moment decreases continually to approximately 
120 degrees angle of attack beyond which the pitching moment 
increases to zero at 180 degrees angle of attack. 

In addition to showing body plus fin data, Figures 
19 and 20 also show typical installed fin data. The measured 
quantities associated with the installed fin data are the 
same as those measured for the fin alone data. The effect of 
the presence of the body on the variation of fin normal force 
with angle of attack is a function of span ratio at angles of 
attack below 75 degrees. For span ratios of 0.3 and 0.4 the 
effect is to slightly decrease the level of fin normal force 
while for a span ratio of 0.5, the level is either maintained 


or increased slightly. For all fins, the presence of the body 
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causes an increase in the level of fin normal force at angles 
of attack above 75 degrees. The effect of the presence of 
the body on fin normal force was consistent throughout the 
Mach number range. 

Another significant effect of the presence of the 
body is seen in Ycpp,,, the lateral center of pressure of 
the fin. Since the spanwise pressure distribution of the 
fin alone is modified by the presence of the body the effect 
ig to move the center of pressure of the installed fin in- 
board over the complete angle of attack range. This trend 
holds true for all of the fins at all Mach numbers. Net only 
does the body alter the pressure distribution over the fin but 
the fin also aiters the pressure distribution over the body. 
This modification of the fin pressure distribution by the body 
will later be referred to as the body on fin interference ind 
the modification of the body pressure distribution by the fin 


will be referred to as the fin on body interference. 


Support Interference 

The mismatch in the body alone normal force at 
90 degrees angle of attack and M = 0.8, Figure 23, for the 
nose mounted and aft mounted struts was an indication of 
possible support interference. Another indication was the 
disagreement between data obtained using a sting support and 
data obtained using a strut support at angles of attack 


between 70 and 90 degrees, Reference (1). In order to 
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determine whether or not support interference was present, 

a series of tests was conducted including free flight aero- 
ballistics range tests at transonic speeds and angles of 

attack of 90 degrees. It was determined that. the aft strut 
supported configurations gave measured normal-force coefficients 
which were too low at Mach numbers of 0.6 and 0.8 and angles 

of attack from 70 to 90 degrees, with litte or no support 
interference indicated at higher Mach numbers. A second 


series of tests, Altstatt and Dietz (51), used sting supported 


models with dummy struts and strut supported models with 
dummy stings to determine the extent and the magnitude of the 
support interference for the 2£/d = 10 body alone configuration 
from the AEDC high angle of attack data set Reference (1). 

The normal-force coefficient corrected for support inter- 
ference is shown in Figure 23 for the 2£/d = 10 configuration. 
Only the data over the angle of attack range from 70 to 90 
degrees, obtained with aft mounted strut model of Figure 8a 


was corrected. The data obtained using the nose mounted strut 
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model of Figure 8b over the angle of attack range from 90 
to 180 degrees matches exactly with the corrected data, 
indicating that there is little or no support interference 
associated with the nose mounted strut. The nose mounted 
strut does not interfer with the body wake whereas the aft 
mounted strut is located on the oody and lies in the body 


wake acting ike a splitter plate, reducing the crossflow 
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drag of the body and thus reducing the normal force, 
2 Nelson (52). 
For the configurations longer than £/d = 10 in | : 
Figure 23, the effect of the strut cannot be determined : 
exactly; however, it would be expected that some degree 
of support interference exists over the same angle of 
i attack range. For the lorger models, the support strut { 
: ' intersects the body in both the aft and forward mounting | | 
| configurations used to obtain data in the angle of attack ‘ 
: range from 0 to 90 degrees and 90 to 180 degrees, i 
t respectively. For these configurations, the strut is always i 
‘ in the body wake. : 
: There was no noticeable effect of the support on the ; 
measured pitching-moment coefficient at the angles of attack ; j 
and Mach numbers at which the tests were conducted. 
: in contrast with the body alone data of Figure 23 i 
where support interference was indisated by the mismatch i 
: of data at 90 degrees angle of attack for the two support | 
; strut configurations, the body plus fin data of Figure 25 i 
. show no mismatch at 90 degrees. In Figure 25, the nose mounted 

: strut data blend smoothly with the data from the aft-mounted 

strut. Since the addition of fins to the model for both 

the nose mounted and aft mounted strut configurations provide 

the same measured normal force at 90 degrees angle of attack, 

it is possible that the presence of the fin on the leeside 

: } 


of the model has an effect on the body wake similar to that 
of the strut. If this were the case, then the effect of the 
strut would be reduced or eliminated. A better understanding 
of support interference for finned slender bodies at high 


angles of attack is needed. 
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V. SEMI-EMPIRICAL THEORY 


The prediction of normal force and. pitching moment 


coefficients for slender body configurations using the 
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Slender body aerodynamics theories discussed in Section II, 
while adequate for preliminary design at low angles of at- : 


‘ tack are generally not adequate at angles of attack above 
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approximately 45 degrees. A semi-empirical theory based on : 
a modified crossflow theory with empirical relations for the 
effects of tail fins, which is adequate for preliminary de- 


sign purposes, is herein proposed. 
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The total normal force and pitching moment for a 


given slender body configuration is made up cf contributions 


wee Pee 


by each component of the configuration and the mutual inter- 


ference of the components with each other. The dominant. con- 
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tribution comes from the body alone forces and momenta with 
the fin alone ferces and the fin center of pressure reitative 

‘ to the moment reference center providing the next largest con- 
tribution. ‘The interference of the body on the fin normal 
force and the interference of the fin on the body normal 

force contribute to the total normal force, while the inter- 
ference forces along with their effective centers of pressure 
contribute to the total pitching moment. The interference 


forces on the fin and body are those which result from the 
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modification of the pressure distribution on one component due 
to the presence of the other component. The determination of 
the interference forces and their centers of pressure at. small 
angles of attack is discussed by Pitts, et al. (53). The 
total normal force on the configuration is described in Refer- 


ence (53) as being made up of a linear combination of the 


component forces and the interference forces. The total pitch- 


ing moment on the configuration is described as being made up 
of the component forces acting at their centers of pressure 
and the interference forces acting at their centers of pres- 
sure. For the empirical determination of the interference 
forces used in theory developed herein, incremental inter- 
ference forces ACy,op, and ACNrop are evaluated frem the 
measured forces on the individual components and combinations 
of components. Effective centers of pressure, Xcppop, and 
XcPpoy, are then eva..uced from the measured moments cf the 
individual and combined components. The centers of pressure 
of the interference forces are described as effective centers 
orf pressure because they are not determined from actual pres- 
sure distributions but only inferred from the measured data. 
The semi~empirical calculation procedure which follows pro- 
vides for the determination of cack contributing factor to 
the buildup of the normal force and pitching moment and then 
eombines the factcers linearly to determine the totai normal 


force and pitching moment ccefficients. 
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It is assumed in the development of this prediction 
technique that the total normal force coefficient, for a 
Slender body with four fins arranged in a cruciform plus 
orientation, is made up of a linear combination of the con- 


tributions of each component given by: 


S¢ S¢ 
Cy = Cupa + 2(CNpa) |B] + ACnpop [= 


S¢ 


where each component is converted to a common reference area. 
The pitching moment coefficient is likewise assumed to be a 
linear combination of the force contributions along with 
their centers of pressure or effective centers of pressure. 


The pitching moment equation for the slender body with four 


fins is given by: 
S¢ S¢ 
Cm = Cmpa * 2(Cupa) (Xcppa) (] + (ACNpop) (Xcppop) | 5 


s 
+ 2(ACNgop) (XcPpop) (2) (5.2) 


where XcPpon and Xcpgop are the effective centers of pressure 
of the incremental forces due to interference. The develop~ 


ment of the calculation procedure is now established. 
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Body Alone Method 
This part of the procedure calculates the forces and 


moments for slender finless bodies at angles of attack to 180 
degrees. The method is based on a modification of the cross- 
flow theory formulated by Jorgensen (6). The crossflow drag 
coefficient variation with Mach number and Reynolds number is 
a modification of the variation reported by Fidler and 
Bateman (8) and the variation ueed in the USAF Datcom (26). 
The equation for normal force coefficient for the angle of 


attack range 0 £ a £ 180 degrees is given by: 


Ss - 
Crea = (2) sin (2a°) cos (3°) 
Ss 2 , 
+ nC. (=) sin” (a7) (5.3) 


The modified pitching moment equation for the angle 
of attack range 0 £ a £ 90 degrees is given by: 


Vv ™ S (2-X. ) a 
Cmpa = [=| sin (2a*) cos (=) 


+ nCa, + ——y— sin” (a*) +2 (5.4) 


and the pitching moment equation for the angle of attack 


range 90 < a £ 180 degrees is given by: 


58 


GTI EAT AE ES ORE EE Ee A Fe 


) 


9 Se er eee Ge eee ee See ere oe 


A GA ER IE SETH AA EER, COLT IONE RE ee WLI Oe SE EE nr ee 


LONE GR LEC IRR OIY TS DINE MSRM HIP 


2 ER SES AER Te EEE IL MPS TERR NER TA 


ea eRe The gerne ar 


Vv- S - 
Ca, * 7 lI] sin (20°) cos ($) 


s (X, = X) 
+ nea. (2) [=—] sin? (a°) + 2 


@ 


where a* =a 0 Sa £ 90 degrees 


a° = 180 -a 90 <a < 180 degrees 


The location of the aerodynamic center given by: 


Cmpa 


x So eee 
“CPRA CNBA 


(5.5) 


(5.6) 


(5.7) 


is measured from the moment reference location and is positive 


forward of the moment reference point. 


The term n is used to modify the two-dimensional 


drag coefficient to approximate the drag coefficient for a 


finite length cylinder and is determined from the datu ob- 


tained by Goldstein (7). The variation of n as a function of 


length to diameter ratio of the cylinder, shown in Figure 28, 


was obtained by a least squares, fifth order polynomial curve 


fit approximating Goldstein's n curve in Reference (7). 


It 


has been customary in the past to assume that the finite length 


correction applied only at subsonic Mach numbers and at Mach 


numbers of 1.0 or greater, the term was constant and equal to 


1.0. This assumption causes a discontinuous change in the 


normal force and pitching moment at a Mach number of 1.0. 
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However, since there is a mixture of both subsonic and super- 
sonic flow over the body at high subsonic Mach numbers, a 
rapid but smooth change in n would be the most likely varia- 
tion. Also the discontinuous change in n at M= 1.0 results 
in an overprediction of Cy at M = 1.0 and a = 90 degrees. 
Therefore, a hyperbolic tangent variation in n over the 
region 0.95 < M < 1.35 has been assumed in the development of 
this technique. The variation of n with 2/d of the configura- 
tion is determined from Figure 21 and then modified for Mach 


number in the range 0.95 £M £ 1.35 by the following equation: 


n= n+ [(1.0 = n)/2] [1.0 + Tanh{ (M-1.0) (15.0/M*)}] (5.8) 


The fourth power of the Mach number in the last term allows 
for a rapid increase in n in the region 0.95 £ M$ 1.0 with a 
slower increage in n in the region 1.0 <M S$ 1.35. 

The crossflow drag coefficient, Cagr used in this 
procedure, is shown in Figure 29 as a function of crossflow 
Mach number, free-stream Mach number, and crossflow Reynolds 
number. For a crossflow Mach number up to 0.6, the cross- 
flow drag variation was taken from Reference (8). At higher 
crossflow Mach numbers, the crossflow drag is assumed to be 


a function of crossflow Mach number only and is represented 


by a modification of the crossflow drag from the USAF 
Datcom (26). 
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The term Z, which appears in the pitching moment 
equation, is the empirical modification to the crossflow 
theory to make the theory fit the observed vitching moment, 
data from the high angle of attack data base. The term is 
a function of both Mach number and angle of attack for Mach 
numbers less than 2.0. The term Z was determined by sub- 
tracting the pitching moment coefficient calculated by 
Jorgensen's formulation of the crossflow theory from the 


measured pitching moment coefficient: 


2 = Cineasured ~ “Mealculated (5.9) 
For each Mach number, 2, as a function of angla of attack, 
was normalized by its maximunm'‘value. ‘The resulting 6 is 
shown in Figure 30 for each Mach number. A curve, 3 (weighted 
toward M = 0.9), also shown in Figure 30, represents the varia- 
tion of 6 with angle of attack. The 5 was represented by a 
Chebyshev polynomial for machine computation. The normalizing 
factor, Zmay, is shown in Figure 31 as a function of Mach number. 
The Zmax variation with Mach number was also represented by a 
Chebyshev polynomial for machine computation. The coefficients 


of the polynomial for both & and Zyay are given in Table 5. 


Another empirical factor was added to account for 
a body with (2/d) different from the data used to determine the 


correction. Thus the resulting factor is given by: 


2 
t/a 
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Coefficients of Chebyshev Polynomials 
‘ Fee na Ee EE TD 4 
z 3 MAX 
| -1.19012E-01 5 56792E+00 
; -2.61068E-01 ~2,86348E+00 : 
[ 1.41951B-01 -3.83133E-01 i 
; 4.53981E-01 2.15263E+00 
-3.71147E-02 4.52087E-01 : 
f ~2.38234E-01 ~3.49564E-01 
8.02296E-03 3.80121E-01 i 
2.41076E-02 ~1.24466E-01 
-2.44889E-02 ~1.07408E+00 
2.18215E-02 -6 .48346E-C2 
5.000118-03 ~1.05244E+00 j 
-3.46462E~-02 ~3.29285E-01 1 
2.56479E-02 ~9,00879E-01 
2.18353E-02 F 
-4.902955-01 
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j 
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Interference Factors 


The incremental interference force coefficients 
ACNrOR and ACNRor: or interference factors, along with their 
effective centers of pressure were determined from the data 
in the high angle of attack data set. Only the data obtained 
for the t/d = 10 total length configurations were used to 
empirically determine the interference factors. The computer 
program used to determine the interference coefficients is 
described in Appendix D. 

The incremental normal force coefficient on the fin 
due to the presence of the body was determined by subtracting 
the normal force coefficient measured on the fins alone from 
the normal force coefficient measured on the fins in the 


presence of the body. 


The above interference factor has a reference area based on the 
fin area and will be converted to a reference area based on 
body cross sectional area in the final normal force and pitch- 
ing moment equations. 

- The incremental interference normal force coefficient 
on the body due to the presence of the fin can now be deter- 
mined by rearranging the assumed equation for the total 


normal force coefficient, Equation 5.1. 
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Again this interference is based on fin area and will be con- 
verted in the final equation. Now with the two interference 
factors known, their effective centers of pressure must be 
determined. The effective center of pressure in the X direc- 
tion relative to the fin hingeline can be determined for the 
interference factor ACNROF by using the measured fin hinge 
moment fer the fin alone and fin in the presence of the body 
cases. The interference hinge moment is obtained by sub- 
tracting the hinge moment measured on the fins alone from the 
hinge moment measured on the fins in the presence of the 
body. The following equation is given for the interference 


hinge moment 
resulting in 


(5.4) 
NBOF 


XCPary ™ 


The lateral effective center of pressure YCPaor can be deter- 


mined in a like manner from the measured fin alone and 
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installed fin root bending moments. 

The effective center of pressure of the interference 
is related to the hingeline and is nondimensionalized by the 
root chord length. It must now be determined relative to the 
center of gravity of the configuration and nondimensionalized 
by the body diameter for use in the final pitching moment 
equation. The hingeline location, Xun? relative to the center 
of gravity for finned bodies, is an input parameter and is 
negative aft of the center of gravity. It follows that: 


Cc 
Xcppor * %w, * Xcrarn [4 (5.15) 


The effective center of pressure of the interfer- 
ence factor, 4CNrop’ can now be determined by rearranging the 
equation for the total pitching moment coefficient, Equation 
5.2. 


s s 
Co Cmpa” 2 (Cre) [+] cesar) ~ 2(CNpa) XcPgq) 3| 
SCirop | 3 


(5.16) 


*cProg 
where 


Cc 
Xcpra "= Xu + CPxna (z] (5.17) 


By the nature of the equation, the effective center of pres- 
sure, Xcprog, i8 nondimensionalized by the body diameter. 
The interference factors and their effective centers 


of pressure have been mathematically represented by a 


65 


hypersurface determined at each Mach number and angle of 
attack, using a multiple linear reqression technique. The 
dependent variable on the surface ja represented as a functicn 
of the three ratios - taper, mapack, sna Span - which define 


the fin. 


Fin Alone Method 

The fin alone contribution to the total normal force 
coefficient and total: pitching moment. coefficient is deter- 
mined by a surface fit to the measured data at each angle of 
attack ard Mach number. The surface was determined ina 
manner similar to the interference factors by the multiple 
linear regression technique. The fin alone variabies of fin = 
normal force, Cigar center of pressure location in the xX 
direction, CPy,,,, and center of pressure location :n the Y 
direction, CPypoqr are determined at. each Mach number and 
angle of attack by an equation which is a functioh of the 


two ratios, taper and aspect, which are independent of the 


body and define the fin. 


Multiple Linear Regression Technique 

The multiple linear regression technique, used to 
represent the calculated interference factors and the 
measured fin alone data by surface equations, is a standard 
application program for the IBM scientific subroutine packaye, 


Yeference (54). 
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Subroutines from the scientific package used to 
perform the multiple linear regression are CORRE, ORDER, MINV, 
and MULTR and a detailed description of each may be found in 
Reference (54). These subroutines were incorporated into an 
overall program described in Appendix E, which prepared the 
calculated coefficients for analysis by defining the surface 
equations and setting up matrices containing the dependent 
and independent variables. A regression analysis was per- 
formed for each parameter at each Mach number and each angle 
of attack combination. Since ten Mach nunhers and 91 angles 
of attack were used, 910 regression analyses were conducted 
for each of the eight parameters, AOCNeop, 4CNpop- XcProp: 
YcPpor’ XcPgry’ Nay CPxXyLar and CPypca: 

Subroutine CORRE was used to determine the means, 
the standard deviations, and the correlation matrix for the 
parameter being analyzed. The subroutine ORDER then selected 
a dependent variable and a svbset cf independent variables 
from the larger set of variables resulting from subroutine 
CORRE for the parameter. ‘The correlation matrix of the subset 
selected by ORDER was inverted by subroutine MINV and finally 
the regression coefficients and confidence level indicators 
were determined by subroutine MULTR. The regression coeffi- 
cients for each parameter and the equation for that parameter 
are takuiated in Appendix B tor the interference coefficients 


and Appendix C for the fin alone parameters. 
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vI. DISCUSSION OF COMPUTATION TECHNIQUE 
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The calculation of the normal force and pitching 


palin au wa 


monent for a finned slender body may be carried out by either 
hané or ¢achine computation. An example of each will be pro- 
vided. Ths hand eaiculation will provide estimations where 
only a few angler of attack and Mach raunmber cases are needed. 
Where more extencive calculations are required, the computer 
program will provide yguick answers for a minimum amount of in- 
put. The range of inputs to the computation technique is 


shown ir. Table 6. 
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Hand Calculatica 

For a given finned slender body, the total normal 
force and pitching moment coefficients are given by Equations 
5.1 and 5.2 respectively. Each component of the two equations 
is determined separately. The body alone normal force is 


determined from Equation 5.3 and the body alone pitching mom- 
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ant is determined by Equations =.4 and 5.5. The factors Cag 
and n are determined from Figures 29 and 28 respectively. The 
factor Z is determined from Equation 5.10 with ¢ and 2 AX 
: determined from Figures 30 and 31, respectively. 
The interference factors and their effective centers 
of pressure are determined from the regression coefficients, 


Bo-++B,g, tabulated in Appendix B, and the general equation: 
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Factcr = By + B,A + Bar? + B, AR + 8, id/b') (6.1) 
The fin alone centributions are determined by the regression 


coefficients, 8)...B,, tabulated in Appendix C, and the 


general equation: 


Fin Alone ~ 


Ry + 6 


2 
3» + BoA + By AR (6.2) 


Since the regression coefficionts are determined for discrete 


Mach numbera and angles of attack, values for other Mach num- 


bers and angles of attack mist be determined by linear inter- 
polation. An example of a hand calculation is given in 


Appendix F. 


Machine Calculation 

Machine computations are carried out using an IBM 
370/165 computer. The program used is described in Appendix 
Ge Only the salient features of the program will be de- 
scribed here. The basic descrirtion cf the nose, body, and 
fins are input as well as the fliyht conditions such as Mach 
number and either altitude or Reynolds number. A subroutine 
in the program calculates Reynolds number from Mach number 
and altitude using data from Reference (55) if altitude rather 
than Reynolds number 18 input. The equations from Appendix I 
are programmed so that only the basic dimensions of nose type, 
nose length, nose radius, body diamecer, and body length are 


required to establish the slender hody configuration. The 
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span ratio, aspect ratio, and taper ratio describe the fin 
configuration. 

The reyression coefficients are input to the program 
from magnetic tape and are called into the program and stored 
on a disk file for use when needed. For each configuration 
input to the program, the aerodynamic coefficients are cal- 
culated for each of the ten Mach numbers: 0.6, 0.8, 0.9, 1.0, 
1.15, 1.3, 1.5, 2.0, 2.5, and 3.0. In addition, for each Mach 
number, the coefficients are determined at angles of attack 
from 0 to 180 degrees at 2-degree intervals. After all of the 
coefficients are calculated, linear interpolations are per- 
formed at each angle of attack to provide calculations at de- 
sired Mach numbers other than those used in the primary calcu- 
lations. For each desired Mach number, the coefficients are 
printed for each angle cf attack from 0 to 180 degrees at 
2-degree intervals. 

The effects on the normal force and pitching moment 
coefficients of roll of the fins from the vertical and 
horizontal planes to an arbitrary roll position, ¢ is accom- 
plished to a first order approximation by multiplying the 
terms containing Cura and ACNE oF in the normal force and 
pitching moment coefficient, Equations 5.1 and 5.2 respec- 
tively by (sin ¢ + cos $). It should be noted that the 


effects of the fins rolling through the body vortices is not 


included in the approximation. The fin normal force coefficient, 
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root bending and hinge moment are determined for the ¢ * 0 i 3 
| case only. 
| Machine caiculations were used for the calculation 


of coefficients for comparison with data for the verification 
of the computation technique. An example of a machine cal- 


| 
| culation is given in Appendix H. 
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VII. VERIFICATION 


Due to the almost total lack of unclassified high 
angle of attack data for finned slender bodies at transonic 
Mach numbers, the verification of the computation technique 
will have to be based primarily on data from the high angle 
of attack data base. The verification will be accomplished 
in three ways. First, computations made using the aerodynamic 
coefficient prediction technique developed herein, will be 
compared with typical data from the high angle of attack data 
base which were used to determine the interference factors. 
Second, computations will be compared with data from the high 
angle of attack data base, which were not used to determine 
the interference factors. Finally, computations will be com-~ 
pared with some of the limited amount of unclassified, high 


angle of attack data from the literature. 


High Angle of Attack Data Used to Determine Interference 
Wetrictents 


Comparisons are made in Figures 32, 33, and 34 
at Mach numbers of 0.8, 1.3, and 2.0, respectively, of the 
measured and calculated aerodynamic force and moment coef- 
ficients for the £/d = 10 ogive cylinder body used as the 
body alone configuratic: for determining the interference 
coefficients. While data are compared only at one 


subsonic, one low supersonic, and one supersonic Mach 
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number, they are typicai of the comparisons at tha other Mach 
numbers. The mismatch in the measured and calculated center 
of pressure at very low and very high angles of attack is the 
result of missing the pitching moment slightly and dividing 
by a very small normal force coefficient. Since the symbols 
on the normal force and pitching moment plots approximate the 
error band of the data, it can be seen that at M = 0.8, the 
calculated values of normal force coefficient lie within the 
error band for most of the angle of attack range, except 
around 90 degrees, where the modification of Jorgensen's body 
alone crossflow theory does not account for the reduced 
normal force when the wake becomes unsteady and where the 
support interference has affected the data. The pitching 
moment coefficient prediction is within the error band of the 
data for much of the angle of attack range except for two 
segments around 55 degrees and from 135 degrees to 165 de- 
grees. At M = 1.3, thenormal force coefficient prediction is 
within the error band of the data except for a small segment 
around 90 degrees that is within a maximum of approximately 
15% of the measured data. The pitching noment coefficient 
prediction is within the error band of the data over much of 
the argle of attack range, with a maximum deviation of ap- 
proximately 15% at 155 degrees angle of attack. At M = 2.0, 
the normal force coefficient prediction is within the error 
band of the data over mostof the range except from 100 te 149 


degrees where the maximum deviation is less than 10% of the 
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Measured data. The pitching-moment coefficient prediction 
is on the edge of the error band over most of the angle of 
attack range except from an angle of attack from 109 to 140 
degrees where the maximum deviation is approximately 12% 
from the measured data. 

A typicai finned configuration used to determine 
the interference coefficients is configuration N2A1T35. The 
data for this configuration are compared with predicted 
values at Mach numbers of 0.8, 1.3, and 2.0 in Figures 35, 
36, and 37, respectively. The center of pressure, normal 
force coefficient and pitching moment coefficient agree with- 
in approximately 15% over the entire angle of attack range 
at M= 0.8. For finned configurations, the installed fin 
characteristics such as fin normal force coefficient and the 
X and Y centers of pressure are calculated. Also determined 
from the fin normal force and the centers of pressure are 
the hinge moment and root bending coefficients for the fin. 
The predicted fin normal force coefficient in Figure 35 is 
within the error band of the data except for the peak due to 
vortex-lift on the triangular fin at approximately 30 degrees. 
The pronounced peak does not occur for the A = 1.0 and 0.5 
fins and was not effectively represented in the surface fit 
since the variation with taper ratio was represented by only 
a second order form in the equation. The root bending and 
hinge moment are represented within the error band of the 


data over much of the angle of attack range and fairly well 
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over the rest. of the range. The X and Y centers of pressure 


La eee sake Lael at st 


show the same variation from the data that the fin moments 1 
show. The comparison of the measured and predicted aerody- 
namic coefficients for the configuration N2A1T35 at a low 


supersonic Mach number is shown in Figure 36 at M = 1.3. The 
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center of pressure and pitching moment coefficient predic- 
tions are within the error band of the data over much of the 
angle of attack range, with maximum deviations of approxi- 
mately 12%. The normal force coefficient prediction is with- 


in 10% over the angle of attack range where it is not within 
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the error band of the data. For this case the fin normal 
force coefficient prediction is within the error band of the 
data. The hinge moment prediction is in close agreement 
with the data except between 35 and 45 degrees angle of at- 
tack and again at about 15G degrees. The root bending is 
also in close agreement with the data. The X and Y centers 
of pressure again show thesame variation from the data as 


the hinge moment and root bending. whe discrepancies in the 
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X and Y centers of pressure are the result of a poor surface 
fit of the calculated value of the effective center of pres- 
sure of the incremental interference force cn the fin. These 
discrepancies do not detract from the technique as a predic- 
tive tocl because of the nature of the X and Y center of 
pressure variation with angle of attack and the resulting 
hinge moment and root bending variation with angle of attack. 


Since the discrepancies ovcur in a region of fairly constant 
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variation of the parameters with angle of attack, a line 
could be faired through the mean of the discrepancies and 
result in a prediction within less than 10% of the measured 
Jvalues. The predicted normal force at M = 2.0 lies within 
the error band of the data at angles of attack up to 100 de- 
grees; beyond 100 degrees the maximum deviation is less than 
10% vrom the measured normal force. The pitching-moment co- 
efficient which was predicted lies within the error band ex- 
cept for a region between 90 and 140 degrees where the pre- 
diction follows the trend of the data with less than 14% 
deviation from the measured pitching moment. Excellent agree- 
ment is obtained for the fin variables with the same type de- 


viations discussed for the M = 1.3 case. 


Other Data from High Angle of Attack Data Base 

It would be expected that the data used to develop 
an empirical prediction technique would ke in excellent agree- 
ment with the predicted values as it was in the previous sec- 
tion. In this section, data from the high angle of attack 
data base using the same fins but a different length body 
will be compared with coefficients calculated using the high 
angle of attack coefficient prediction technique. 

The measured aerodynamic coefficients for the t/d = 
15 bedy alone configuration are compared with predicted co- 
efficients in Figure 38 at a Mach number of 0.8. Just as in 


the subsonic data for the t/d = 10 configuration, excellent 
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agreement is obtained except for the ar. le of attack range 
from 60 to 120 degrees where the unsteady wake and suspected 
support interference affect the measured data. Agreement 
within 15% is obtained in the pitching moment coefficient ex- 
cept in the angle of attack range from 120 to 160 degrees 
where the technique overpredicts the pitching moment coeffi- 
cient by a maximum of abcut 40%. The center of pressure is 
in error for angles of attack above 120 degrees due to the 
error in pitching moment. 

Suversonically, as shown in Figure 39, at a Mach 
aumber of 1.3 the predicted aerodynamic coefficients for the 
2/d = 15 body alone configuration are in much closer agree- 
ment with the data. The norma] force coefficient is in 
error by less than 10%. The correction to the pitching mo- 
ment coefficient derived for the 2/d = 10 configuration 
causes a flattening in the curve between about 20 and 50 de- 
grees angle of attack, resulting in the underprediction of 
the pitching moment coefficient in this range. Over the rest 
of the angle of attack range, agreement within a maximum of 
15% is obtained. 

The subsonic case at M = 0.8 for the finned L£/d = 
15 configuration, N2A3T31 with AR = 0.5 fins, is shown in 
Figure 40. The predicted nérmal. force coefficient is within 
a maximum of 15% of the measured value and except for one 
small range of angles about 33 degrees, where an error of ap- 


proximately 20% occurs, the pitching moment coefficient 
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prediction is within 15% of the measured value. Most of the 
error in the coefficients can be attributed to the error in 
the predicted body alone coefficients. The predicted, in- 
stalled fin normal force agrees with the measured values 
within less than 15% over the complete angle of attack range. 
The x and Y¥ centers of pressure and the fin moments again 
display the same discrepancies discussed earlier but again 

A the level of the parameter can be estimated by fairing. 


The supersonic case for configuration N2A3T31 is 
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shown in Figure 41 for a Mach number of 1.3. Excellent agree- 
ment is seen for this configuration between the measured and 
predicted normal force pitching moment and installed fin 
normal force coefficients. The coefficients are predicted 
within 10% over the entire angle of attack range. Even the 
predicted fin moments and centers of pressure are in close 
agreement with the measured values. 

The data for a second finned 2/d = 15 configura- 
tion, N2A3T13, with AR = 2.0 fins, is now compared with the 
predicted coefficients. The subsonic case for M = 0.8 is 
shown in Figure 42. As in cther cases at subsonic Mach num- 
bers, excellent agreement in normal force coefficient is ob- 


tained except in the unsteady wake region around 90 degrees 
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angle of attack where a maximum error of approximately 15% 
occurs. The pitching moment coefficient agrees almost with- 


in the error band of the data except at angles of attack 
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between 116 and 146 degrees where a maximum error of approxi- 
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nately 206% occurs. The fin characteristics of normal force 
and root bending are slightly underpredicted at angles of at- 
tack above 80 degrees while the other fin characteristics are 
in excellent agreement. Complete data were not available at 
M= 1.3, out the comparison of measured and predicted covet ‘i- 
cients at M = 1.15 is shown in Figure 43. The measurements 
and predictions are in excellent agreement, which is typical 


of the previous comparisons at supersonic Mach numbers. 


High Angle of Attack Data from Literature 


In this section, data from recent tests at high 
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angles of attack are compared with the theoretical predic- 


tions. Data for a slender body repo. ted by Fleeman and 
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Nelson (36) provide a comparisun of the body alone predictive 
capability of the prediction technique. The slender hody 
consists of a 2.5 caliber sharp ogive nose with a 12.05 cali- 
f ker cylinder afterbedy. The model was tested tsing a sting 
Support at angles of attack from 0 degrees to 45 degrees and 


from 135 degrees to 160 degrees. At angles of attack trom 


45 degrees to 135 degrees, a strut support was used. The 
strut support intersected the aft portion of the body ata 
90 degree angle. The predicted vaiues of normal force and 


pitching moment are ccmpared with the measured daca in Fig- 
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ure 44. At Mach numbers of 0.6 and 0.8, the normal force 
agreement is excellent except around 90 degrees angle of at- 


tack, where th nsteady wake effect and possible support 
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interference cause approximately 17 to 20% disagreement. 

The predicted pitching moment coefficient agrees with the 
data for both M = 0.6 and 0.8 up to 100 degrees angle cé at- 
tack. Between 100 degrees and 160 degrees, the theory pre- 
dicts a much lower value of the coefficient than was meas~ 
ured. Just as with the previously described data, the agree- 
ment is excellent at the supersonic and low supersonic Mach 
numbers. 

Another slender body configuration test was con- 
ducted by Baker and Reichenau (40). The tests were conducted 
for a series of air and ground launched strategic missile 
concepts. The data for a typical configuration from the test 
are compared with predicted normal force and pitching moment 
coefficients in Figure 45. The configuration selected is 
N3B2. The N3 designates a blunt ogive nose, 2.14 calibers in 
length with spherical blunting. The afterbody, B2, is a 


cylinder 6.15 calibers in tergth. The model was supported 


by .a sting in the angle of attack range from 0 degrees to 45 


degrezs and with a forward swept strut in the angle. at-~ 
tack range from 40 degrees to 180 degrees. The large forward 
swept strut would be expected to produce more interference to 
the data than the strut of Reference (36), which intersected 
the body at 90 degrees or the support arrangement used in 

the high angle of attack data base which has a 45 degree 
rearward sweep. The discontinuities in the data resulting 


from support interference at 40 degrees angle of attack are 
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evident in Figure 45, where the measured and predicted values 
of normal force and pitching moment coefficient are compared. 
At a Mach number of 0.6, only fair agreement is obtained with 
the normal force coefficient, while very good agreement is ob- 
tained in pitching moment coefficient. Much better agreement 
is obtained at Mach number 0.8 with the error in normal force 
of only about 10%. The agreement in pitching moment is very 
good with greatest disagreement being in the angle of attack 
range from 140 degrees to 160 degrees. At supersonic Mach 
numbers the effect of the strut appears to cause a decrease 
in the level of the normal force and pitching moment coeffi- 
cient. The agreement between the measured and predicted 
normal force and pitching moment coefficients is excellent 
for the data obtained with the sting support; however, the 
data from strut support model is in disagreement with the 
theory. 

The only high angle of attack data for a slender 
hody with low aspect ratio fins at transonic Mach numbers 
other than the data in the high angle of attack data base are 
reported by Jenke (47). The tests were conducted on a modi- 
fied basic finner model to measure the roll damping, the 
Magnus force and the static stability. The model was sup- 
ported by an "L" shaped strut arrangement. A six~compenent 
balance was attached to the horizontal leg of the strut and 
a roll mechanism was attached to the balance. The model was 


then attached to the roll mechanism. With the model in place, 


82 


fecha SR aowthe hs 


ne Sl inde ltt able me 


j 
4 
4 
' 
| 
{ 
: 
q 


a a a el a la 


3 


i 
i 


poe ee re rg Tne Soe re Serre SRT nee I EE OT REIS ARNE 


the vertical leg of tne strut was approximately 3.25 medel 
diameters behind the model. The model consis<ad of a 2.5. 
caliber sharp ogive nose an? a 7.5 caliber cylindez aSterbody. 
Four fins were monnted ina cruciform plus orientation. ‘Ine 
fins hac a taper ratio of 0.5, an asnect ratio of 1.0, and a 
gpan ratio of 0.5. The model had a rell cate of approxi- 
mately i100 zsadians/reecond 4uring the test. 

The comparison of the measvred and predicted aero- 
dynamic vcafficients is shown iu Figure 45. At a Mach num- 
ber of 0.6, the maximum disagreement in normal force coeffi-~ 
cient of approximately 36% occurs at 45 degrees. Better 
agreement is obtained at lower angles and at 90 degrees. The 
trend of the pitching moment is correct, including the in- 
crease in the coefficient to a maxinum at approximately 60 to 
65 degrees angle of attack. However, the magnitude of the 
pitching moment at the maximum is not in agreement. The 
agreement in normal force coefficient at a Mach number of 0.8 
is much better than at 0.6, with the predisted values within 
less than 16% of the measured data. Again, the magnitude of 
the measured and predicted walues of pitching moment coeffi- 
cient at the maximum between 30 and 70 degrees angle of at- 
tack do not agree. Supersonically, at Mach numbers of 1.15, 
1.3, 1.5, 2.0 and 2.5, the agreement in both nermal force and 


pitching moment is very good. 
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VIII. CONCLUSIONS 


Conclusions 

A series of slender body configurations, both with 
an@ without tail fins, were tested at transonic and supersonic 
Mach numbers and at angles of attack to 180 degrees. Addi- 
tionally, the tail fins were tested alone on a reflection 
Piane at transonic and supersonic Mach numbers and at angles 
of attack to 180 degrees. The parametric variation of model 
dength and the fin aspect, taper, and span ratios provide a 
significant data base for body alone, fin alone, body plus 
fin and installed fin configurations. Prior to this effort, 
there was a complete dearth of parametric data for slender 
bodies with low aspect ratio fins at transonic and supersonic 
Mach numbers in the very high angle of attack range. The data 
base established through this effort provides a hasis for the 
development. of the semi-empirical aerodynamic coefficient pre- 
diction technique reported herein. The data base also wiil 
provide a standara of comparison for completely theoretical 
aerodynamic ccefficient prediction techniques being developed. 

A semi-empirical theory, adequate for preliminary 
design purposes, has been developed for the prediction of aero- 
dynamic coefficients for finned slender bodies at angles of 


attack from 0 to 180 degreee and Mach numbers from 0.6 to 3.0. 
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The theory is based on & modification of the crossflow theory 
as formulated by Jorgensen with empirical relations for the 
effects of tail fins. An empirical set of interference fac- 
tors was determined to correct the fin alone data for the pre- 
sence of the body and allow for the determination of the in- 
stalled fin aerodynamic characteristics. Empirical inter- 
ference factors were also determined to account for the incre- 
ment in body force due to the presence of the fin. A multiple 
linear regreesion technique was used to put the vast amount of 
interference factor and fin alone data into a form which is 
simple to use and provides the capability for the determina- 
tion of interference factors and fin alone data for arbitrary 
low aspect ratio fins within the range of the data base. 
Equations involving regression coefficients and the aspect, 
taper, and span ratios of the fins are used to calculate the 
interference factors and fin alone aerodynamic coefficients. 

A computer program has been written to provide 
rapid calculation of the aerodynamic coefficients for multiple 
configurations at user selected Mach numbers and Reynolds 


numbers or altitudes. The aerodynamic coefficients predicted 
by the program are within 15% or better of the measured data 
over most of the angle ot attack range. 
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Figure 25. Typical body plus fin data for three finned 
bodies having APR = 2.0 and d/b' = 0.4 fins with 
different taper ratios, M = 0.8. 
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: Figure 26. Typical body plus fin data for three finned 
bodies having AR = 2.0 and a@/b' = 0.4 fins with 
different taper ratios, M = 1.3. 
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Figuze 27. Tyrical body plus fin data for three finned 
bodies having AR = 2.0 and d/l' = 0.4 fins with 
different taper satios, M= 2.0. 
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hody alone aerodynamic coefficients, M = 1.3. 
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Figure 34. Comparison of typical, measured body alone data 
useé to determine interference coefficients, M= 1.3. 
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coefficients, M = 0.8. 
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APPENDIX D 


DATA ANALYSIS PROGRAM 


The data analysis program 1s relatively short, con- 


sisting of a MAIN program and five subroutines (ORDER, SWAP, 


INTERP, DIM and FORIT). All of the data in the high alpha 


data base, obtained during numerous wind tunnel entries over 
a two-year period were placed on one reel of magnetic tape 


for analysis by the data analysis program. The various data 


part numbers for a given configuration and Mach number which 
go together to make up an angle of attack sweep from 0 to 
180 degrees were given common identification numbers on the 
tape. 

The MAIN program reads the data tape and places the 
data in a sequential disk file to be searched repeatedly to 
bring the various body alone, fin alone and body plus fin 
data into the program so that the interference factors and 
their effective centers of pressure can be calculated. To 
begin the computation for a given fin configuration and a 
given Mach number, the index numbers of the fin alone, the 
body alone and the body plus fin data are read in MAIN. 


The 


disk file 1s then searched for the required data. When one 


of the input index numbers 1s found, the data are placed into 


one of three sets of arrays depending on whether the data are 


aia a 
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ee 
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for fin alone, body alone, or body plus fin. The sets of 


arrays are made up of data for the normal force and pitching 


moment coefficients for both the body alone and body plus fin 


as well as the fin normal force, root bending and hinge mom- 
ent coefficients for both the fin alone and the fin the 


presence of the body. After the search of the disk file is 


complete, the data from each array are placed in an ascending 


order of angle of attack and duplicate angles of attack are 
eliminated by subroutines ORDER and SWAP. Once the arrays 
are ordered, subroutine INTERP 1s called for each array and 


a relaxed second derivative, cubic spline interpolation is 


performed to determine the values of each array at even angles 


of attack from 0 to 180 degrees. The use of the interpolation 


also fills in any small gap in the data due to bad data 
points which were left out or places where problems in the 
wind tunnel test precluded obtaining the data. 

With the data in each array determined at even 
angles of attack, the interference factors and effective 
centers of pressure were obtained using the dimensions of 
the selected fin obtained from subroutine DIM and Equations 
5.11, 5.12, 5.14, and 5.16 from Section V and an array based 
on angle of attack for each factor was established. Because 
of scatter in the data, 1t was necessary to smooth the 
determined values of the interference factors and effective 


centers of pressure as well as the fin alone coefficients 
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before the regression coefficients could be determined. In 
order to do this, the values of sach array were approximated 
by a twenty term Fourier series using subroutine FORIT. Sub- 
routine FORIT is a standard IBM scientific package subroutine 
described in Reference (54). The standard deviation, o, of 
the values in the arrays from the Fourier series approximation 
was then determined. The values of the factors in the arrays 
were then compared with the approximated values using the 
Fourier series. Any value that daviated from the approxime- 
tion by more than 20 was placed equal to the approximated 
value at that angle of attack plus or minus one o depending 
on whether the value from the array was greater or less than 
the approximated value. 

The sinoothed values of the interference factors, 
their effective centers of pressure and the fin alone 
coefficients wera then placed into a larger set of arrays 
which were functionz of angle of attack, Macn number, taper 
ratic, aspect ratio and span ratio. After the values were 


placed into these large arreys control was shifted back to the 


i read statements in MAIN and a new fin configuration was 
analyzed. After all of the fin configurations in the data 


base that were tested with the 2£/d = 10 bcdy had been analyzed 


and placed in the large arrays, the large arrays were then 
written on a magnetic tape for further analysis by the re- 


gression coefficient projram. 
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APPENDIX E 
REGRESSION COEFFICIENT PROGRAM 


The data in the high angle of attack data base 
associated with the fins were obtained for twelve different 
fins each having different combinations of taper, aspect and 
span ratios. A schematic in Figure (E-1) shows the relation- 
ship of each fin to the others in the data base. The three 
ratios describe a cube where one dimension is represented by 
the taper ratio, another dimension is represented by the 
aspect ratio with the final dimension being represented by 


the span ratio. Every combination of the three ratios was 


not tested but the ones where data were obtained are indicated. 


It can be seen that at each angle of attack and Mach number 
combination, a “cube” is defined for each interference factor 
(ACNFog and ACNgop) and effective center of pressure (Xcprop: 
YCPpor’ and Xoppry) - Thus, in order to provide for the 
Getermination of interference factors and effective centers 
of pressure for arbitrary fins with characteristics which 
fall within the "cube", a regression analysis using a least 
squares technique was used to fit a hypersurface to the cal- 
culated factors. The form of the equation for the hyper- 
Surface was assumed from the way the data ir the data base 


were obtained. Since data were obtained at only two 


oe I TT 


Case 


different aspect ratios while holding the other ratios con- 
stant, a linear function of aspect ratio was assumed. Also 
since data were obtained at only two span ratios while hold- 
ing the other ratios constant, t:hen a linear relation was 
assumed for span ratio. For each combination of span and 
aspect ratio, data were obtained for three taper ratios, thus 
a guadratic function of taper ratio was assumed. The hyper- 


surface equation for a typical interference factor therefcre 


is; 


ACNrop = 85 + B,) + 82d” + By AR + By (d/b') 
For every combination of the three ratios at which data were 
taken, an equation is formed. Thus for tne twelve combinations, 
a system of twelve equations is developed. The multiple 
linear regression technique then takes the twelve equations 
and determines the regression coefficients Bors By which 
allows the calculation of ACNrop for arbitrary values of the 
taper, aspect and span ratios. Since the fin alone character- 
istics for a given Mach number, angle of attack combination 
are a function of only taper ratio and aspect ratio, a 


system of nine equations are determined with the form: 


Cup, = By + 8, + BoA? + B AR 


and the regression analysis is performed to determine the 


regression coefficients By: -+B3- 
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The computer program which was written to determine 


peste ne Rae Rie 


the regression coefficients consists of a MAIN program and 
eight subroutines, Figure (E-2), the MAIN program calls sub- 
routine CALCO to prepare the data for analysis. A regression 
analysis is done for each factor at each Mach number, angle 
of attack combination. Thus CALCO sets up a column matrix 
of dependent variables and a 4 x 12 rectangular matrix of 


independent variables which represent the systems of 12 


dich ie ald i oe See us A, 


equations for each of the interference factors and effective 


centers of pressure. CALCO also sets up a column matrix of 
dependent variables and a 3 x 9 rectangular matrix of in- 
dependent variables which represent the systems of nine 
equations for each of the fin alone characteristics. The 
data which make up these matrices are provided by subroutine 
DAP which reads the magnetic tape produced by the data 
analysis program (Appendix D). Since the data are supplied 
by Subroutine DAP, subroutine DATA called by subroutine CORRE 


consists only of return and end statements. After the matrices 


| 
! 
i 
{ 
| 
{ 
| 
| 
{ 
; 
| 
} 


are established for each factor at the first Mach number, 
angle of attack combination subroutine MULREG is called for 
each interference factor, effective center of pressure and 
fin alone characteristic. 

Subroutine MULREG places the column matrix of 
dependent variables into the first column of a rectangular 


working matrix. Then the matrix of independent variables is 
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placed into the remaining columns. With the rectangular 
working matrix subroutine MULREG then calls consecutively 
subroutines CORRE, ORDER, MINV, and MULTR, which were taken 
from the IBM Scientific Subroutine Package, Reference (54), 
to perform the multiple linear regression analysis. Control 
is then returned to subroutine CALCO. 

After subroutine MULREG is called for each factor 
and a regression analysis is performed for each factor, the 
resulting regression coefficients are placed into a large 
array set up for each factor which is a function of Mach num- 


ber and angle of attack. Now a new regression analysis is 


performed for each factor at the next Mach number angle of 
attack combination and these regression coefficients are 
added to the large arrays. When regression analyses are com- 
pleted for each Mach nunber, angle of attack cumbination, 
eight large arrays of regression ccefficients result which 
represent ACNFoR! ACNRoF’ XCProp? YCPpor: and XCPpry 48 
functions of Mach number, angle of attack, taper ratio, 


aspect ratio, and span ratio and CNpae CPxurae and CPyRcA as 


TT Re se OP Te Oe ee eee ee, eee ere Lee ee 


functions of Mach number, angle of attack, taper ratio, and 
aspect ratio. These large arrays are then written on magnetic 


tape to be used in the coefficient pradiction program, and 
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the regression coefficients are printed wut for a permanent 


ee a ON eee ee pe ee 
tae i te ak 


record. 
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Figure E-l. Range of fin variables. 
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Fiaure E-2. Subroutine flow chart for the regression 
ewefficient program. 
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APPENDIX F 
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APPENDIX F 
HAND CALCULATION EXAMPLE 


Calculate the total normal force and pitching 


moment coefficients at a Mach number of 0.6, Reynolds number 

{ 
based on body diameter of 4 x 10° and angle of 120 degrees j 
for a finned slender body having a sharp tangent ogive nose 


and a cylindrical afterbody with the following dimensions: 


Moment reference = 2/2 


(2,/d) = 2.5 
(£,,/4) = 7.5 
a = 1.25 inches 
Xun * -5.5635 inches 


wre noe er 


The four fins are arranged in cruciform plus orientation 


with the trailing edge flush with the aft end of the body 


and have the following characteristics: 


d = 1.0 
AR = 1.0 
a/b' = 0.5 

HL = 0.45 


By hand calculation, Equation 5.4 and the formulas from 


Appendix I determine: 
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Total length £ # 12.5 in. 


Moment ref location X, = 6.25 in. 
Fin area S_ = 0.7813 in. 2 
Base area, Sp = cross 

sectional area, S S = 1.2272 in.? 
Nose-body planform area Sp = 14.3436 in.? 
Nose-body volume V = 13.5738 in.2 
Location of centroid of area X = 6.7391 in. 
Root chord Cr = 1.25 in. 


From Figure 28 for an £/d = 10 configuration, n is determined 


to be: 


n = 0.682 


Since M is less than 0.95, n will not be changed by Equation 


5.8. 
From Figure 29 for M = 0.6 and Re, = 4.17 x 10° 
M. =M sina = 0.5196 
5 
Re, = Rea sin «a = 3.6 x 10 
Cq #2 1.17 
Cc 


Now from Equation 5.3 we get 
From Figure 30 at a = 120 degrees, we get 


ov = -0.88 


and from Figure 31 at M = 0.6, we get 


Thus from Equation 5.10 we get 
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| 
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2 


z= (3) (2y.¥) (443) 
Z = -6.248 

Now from Equation 5.5 we get 
Cma, ™ -11.88 


Now from Appendix B we get the regression coefficients for 
the interference factors and their effective centers of 


pressure. 
ACNpon = 8 * 84) + Bi? + By AR + 8, (d/b') 


where Be m 0.3771 
8, = -0.5460 By = ~-0.0169 
B = 0.5348 By = 0.4997 


ACipor = 8 + 8A + Bah” + By AR + By (A/b’) 
ACNnop = 09-2806 


Where 8B. = 0.3099 


B5 = -0.6732 
By _ -0 e 0338 


B, = 0.2691 
8B. = -0.0567 


Xcppop = 85 + 8,4 + Bg” + By AR + By (d/b') 


XCProp = -2.8492 
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Where 8B. = ~14.1688 
= 2.6929 


4 
| 
85 
8, = 16.2029 


8 = 11.1996 
Bo = -10.6743 


hee, 


, 2 ‘ 
XCPpry * 8, + 8,4 + Bad” + Bz AR + By (d/b') 


et ttl Se 


donee, 


Where B. = -0.0282 
B5 = -0.0294 


arated ae eo 


Bg = 0.0297 | 

Bo = 0.0367 

Yopnor ' 85 + 8,’ + BA” + By AR + B, (a/b") 
Yopgop = 0-0334 | 

Where By = 1.5444 
B, = 0.1914 i Wa 

$i 4 

B = -0.1326 Ba cesaeee : 


From Appendix C we get the regression coefficients to cal- 
culate the fin alone contributions to the normal rorce anda 


pitching moment coefficients. 


lead its Came cn heli Rn os” am ek 


2 
CNrA = By + BA + Boh + By AR 


Cypa = 1.025 

q 

Where B, = 0.8845 B = -0.0228 
B, = 0.1095 B, = 0.0438 
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CPxupa * By + 8,4 + BDA? + By AR 


CPXuTA = -0.1006 
Where Bo = -0.1529 B. = 0.0150 
By = 9.014) Bs = 0.0232 


2 
CP¥pca * 8, + 8,4 + BDA” + BS AR 


CPypca = 0.4705 
Where B, = 0.3528 Bo = ~0.1252 
B, = 9.2502 B, = -0.0073 


The calculated parameters along with the dimensional character- 
istics of the configuration are now substituted into the 
following equations for normal force, pitching moment, fin 
normal force in the presence of the body and the fin X and Y 


centers of pressure in the presence of the body. 


(Ss Se Se \ 


Cy = 9-80 
Xo = xX + CP CR 
PRA HL XHLA +r 


Xcpp, 7-5-6631 


E c 
XcPpor * %yr * XcPpry (| 


XcPpor = ~5-7456 
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S¢ S¢ | 


8 
+ 2(8CNaop) (Xce pop! (7s) 


Cc. = -70.56 


4 

4 

i 

4 
sot 
| eee (Cup, | (CPxuta] . (4Cupor} (craps! | 
HLB Nep | 
; i 
i 
| CPx pp = e+ 118 | 
q 
) (CNea | cPypca} + ( CNpor | ( ¥cPror] 
| CPyacR ® ire 
2 j 
? 
| CPypop * 0-376 

a 
Cmp, = 90-4915 | 
| Cing = 70-154 | 
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APPENDIX G 


COEFFICIENT PREDICTION PROGRAM 


The high angle of attack aerodynamic coefficient 
prediction program is written in FORTRAN IV for the IBM 370 
Computer. It consists of a MAIN program and 14 subroutines 
as shown in Figure (G-1). The MAIN program consists of read 
statements for the title, body geometry, and flight conditicns 
such ac Mach number and either altitude or Reynolds number 
based on body diameter. The title is read using a 20A4 format 
and che geometry for each componert nose, body and boattail 
1s res¢] using an I10 format for the indicator of nose, body 
or boattail type and the lengths and moment reference loca- 
tion are read using an F1C.4 format. The flight conditions 
are read using a format statement containing an implied DO 
loop, where the number of Mach numbers, Reynolds numbers or 
altitudes is read using an 110 format and the Mach numbers 
are read using an F10.4 format. The Reynolds numbers are 
read using an E14.4 format. Either Reynolds number or 
altitude is read into the program and a blank card is inserted 
for the one not read in. Using the input geometry, cross- 
sectional areas are calculated and subroutines OGIVE, NOSE, 
BODY and BTAIL are called to calculate the planform areas, 


volumes and centroids of the various body components using 
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the equations of Appendix I. Subroutine OGIVE determines 
the theoretical length of the og: re if only the ogive radius 
and body diameter are given. The total body planform area, 
volume and centroid are then calculated in MAIN. 

With all of the areas, volumes and lengths deter- 
mined, subroutine COEFF is called to begin the computation of 
the aerodynamic coefficients. For a given body, calculations 
for any number of fins can be made. The number of fins is 
read, using an I10 format. At this point, the fin title, fin 
ratios, and fin orient*tion are read into the program. The 
title is read using a 20A4 format and the taper ratio, aspect 
ratio, and span ratio are read in using an F10.4 format. The 
ratio of the hinge line location measured from the leading 
edge of the roct chord to the root chord length and the 
hinge line location measured in calibers from the moment 
reference point are read using an F10.4 format. The hinge 
line location 1s a negative number if it is located aft of 
the moment reference point. The fin orientation measured 
from the vertical is read in degrees using an F10.4 format. 
If a body alone calculation is desired, then the number of 
fins 1S input as zero and only a fin title card is required. 
A form 18 provided, Figure (G-2), for preparing the program 
inputs. Spaces are provided for five different fin con- 
figurations, but any rsmber can be input. With the fin 


characteristics, subroutine FINDIM is called *o calculate 
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the dimensions and areas of the fins. 

Subroutine PAGER is now called to print a title 
page showing all the dimensions of the body and first fin to 
be calculated. Subroutine CO is now called to determine the 
crossflow drag and the finite length body correction, n. At 
this point, if altitude has been input rather than Reynolds 
number, then subroutine RENOLD is called to determine the 
Reynolds number based on body diameter using the Mach number, 
altitude and curve fits of atmospheric data from Reference 
(55). 

With all of the body dimensions the crossflow drag 
and the n evaluated, subroutine NORM is called to determine 
the body alone normal force and pitching moment coefficients 
and the body alone center of pressure. The empirical correc- 
tion to the body alone pitching moment coefficient: is deter- 
mined by calling subroutine CSPS to evaluate the Chepyshev 
polynomial for the correction. Subroutine CSPS is a standard 
IBM subroutine from Reference (54). 

Now to evaluate the interference factors, their 
effective centers of pressure, and the fin alone character- 
istics, subroutine BAKER is called. This subroutine reads 
the regression coefficients from a magnetic tape and sets up 
a direct access file for subsequent calculations. [It is this 
subroutine that contributes the most to the amount of core re- 


quired to run this program. Since all of the calculations in 
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this program are initially done at Mach numbers of 0.6, 0.8, 
0.9, 1.0, 1.15, 1.3, 1.5, 2.6, 2.5, and 3.0 for even angles 
of attack from 0 to 180 degrees, arrays of interference fac- 
tors and fin alone coefficients are determined as a function 
of the specific Mach numbers and angles of attack using the 
regression coefficients and the taper, aspect. and span ratio 
for the fin configuration being calculated. With the arrays 
of interference factors and fin alone coefficients, subroutine 
DOIT is called to determine the total normal force and pitch- 
ing moment coefficients, the center cf pressure for the 
complete configuration and the installed fin normal force, 
root bending and hinge moment. coefficients as well as the 
X and Y centers of pressure of the normal force on the fin. 
Subroutine PAGER is now called again to print the 
calculated coefficients. Up to this point, ali of the 
calculations have been made at internally determined Mach 
numbers. Therefore, before the results can be printed, 
linearly interpolated values of the coefficients are deter- 
mined using subroutine LINE at the user specified Mach num- 
bers. If input Mach number is below M = 0.6 or above M = 3.0, 
a note is printed and the coefficients are extrapolated t2 
the requested Mach number. Coefficients are printed at even 
angles of attack from 0 to 186 degrees. After the coeffici- 
ents for all of the fin configurations for the first body are 


determined, control is returned to COEFF to start the 
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trol is returned to MAIN to read in the next body configura- 


tion, 


A listing of the Fortran portion of the program 


follows. 
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Figure G-2. Program input form for the coefficient prediction 
program. 
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APPENDIX H 


MACHINE CALCULATION EXAMPLE 


The finned slender body configuration used in the : 
machine calculation example to follow is the same one used 
in the hand calculation example. Included in the exampie 
is the program input form filled out for the example config- 
uration, Figure (H-1), the punched input cards for the 
example configuration, Figure (H-2), the title page print- 
out, Figure (H-3), and the calculated coefficient printout, 


Figure (H-4). 
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Figure H-l1. Program input form for the coefficient predic- 
tion program filled out for the example case. 
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APPENDIX I 
EQUATIONS 
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APPENDIX I 


EQUATIONS 


The following equations are used to calculate the 
plan area, volume and centroid of the nose, body and boat- 


tail components which make up a total body alone configura- 


tion. 


Plan Area of Nose 
Sharp Cone (Figure I-1) 


By = typ (4/2) 


Spnerically Biunted Cone (Figure I-1) 


2 


Sherp Ogive (figure I-2) 


2 2 
bam? + (4/2) ] 


Ay * yo sin ay 


g Q 
zs P - NT we -1 N 
Ay = Lyp (4/2) - Ry [C973 eae Td72) 


2nd 7/2 


4(a/2) | a) y} 


NT“ + (d/2) 


Spherically Plunted Ogive (Figure I-2) 


l. Determine plar area, AN, ! for sharp ogivcs 
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a plate “vie cae tee as oll 


2. Determine parameters, R, Ri, Ro, Lor L3, 


wv. B and 80 
[ 2+ (a/2)* 
R = srry 
ee 
L 
R, = R-R 


us 
{ 
pam 
=z 
ty 
= 
Ro 


Pi = (1/2) (L, + R + Rj) 


Sapo oe > 1/2 
E L3) (Py ~ R) (Py au 


z r, | 
es eart | TEs) 


= [m - 8] 


3. Now determine area to be removed from sharp 


ogive, shaded area (Figure I-1). 
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5.4 
Asa = 3 Ry (28) 


4. The plan area of the spherically blunted 


Ogive is thus 
Ay = Ang ~ Ap t+ Asa 


Plan Area of Body 
Cylinder 


A, = Lp (d) 


Plan Area of Boattail 


Truncated Cone 
dar 
Agr = a (d + aT) 
The plan area of the total slender body is thus: 


Sy = Ay + Ag + Agr 
Volume of Nose 
Sharp Cone (Figure I-1) 
(2/2) 2 (Rin) 
Nyy Seg ee 
Spherically Blunted Cone (Figure I-1) 
1. Determine the parameters Los L3, LL, LLo 


(see Figure I-2) 
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2 NT 3 
2 2 41/2 
(Lyn) “ (Ry) 
2 NT 
Ly = & + lt 
(d/2) 
(d/2) 
a ae SESE 
[car2) ao te ] 


LLy = Lo + Ry 


2. Determine volume from: 


door 


2 2 
_ _1(d/2) 2 uam(LL,) ~— Xam(LL,)° + 
Vy = sg" NT 1 NT 1 


: (LL)? 
+7 R [ ampr- cat] ~1| —e— 


2 (Ly)? 2 


+ (he) oa) 


Sharp Ogive (Figure I-2) 


l. Determine the parameters R, R 


2 2 
_ tee + (4/2) 


Ry = R- (d/2) 


2 
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: 2. Determine volume by: 

| 2 R & 

3. [R- 2. -1 {*wr)\_ 

; Vn =7 Loop E 1.0 - ioe sin (4 1/3 
NT T 


Spherically Blunted Ogive (Figure I-z) 
1. Determine the parameters R, Roe Rie Los LL,, 


LLos Rye Yy- Ripe Ryo Rr 


Cr. 2 
na | + {d/2) 


Pz 
a 1] 
nO 
i) t 

x 6 

™~ 
N 
— 


tc 
nN 
{ 
| een | 
v2) 
Dh 
t 
wn 
Le) 


LL» saad Lo tT Ry 


R = Pee y = Ro R. = oy2 
o£ 1 L LL ~ £0 
1 1 1 
R L 
N _ iL 
Ryo EL; Rh2 ED, 


2. Determine volume by: 


4 
Cae eae ma nea ie ta A ee ct ee te TAR ah i etm a i ANE LOLS Fs Sl a ed ae Man SN an ne a a A uid 
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" 3 2 -1 


N 


2 (R2 ~ 1,0)1/2 


1 
3 


(Ry) 


ah 


64e2 er 
t (Rig — Rpg) (Ry 7 1.0) - SS 


2 
+ Ri (Roy - 1.09] 
Volume of Body 
Cylinder 


Vv. = m(a/2)* (8 


B ) 


B 


Volume of Boattail 


Truncated Cone 


Ver = 5 [ tav2y? + cana)? | (2.,.) 


BT 
The volume of the body is thus: 


VV tv + Ve 


Centroi1d of Nose 


Sharp Cone (Figure I-1,; 

nr, Sr oe 

Xy (273) (kan) 

Spherically Blunted Cone (Figure I-1) 


1. Determine parameter Rin 
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2. Using parameters determined for volume of 


blunt cone, calculate: 


3 
LL 
_ Uy ff (a/2) : 
Ny * Thy am (3 Ring 2.0) 


: to * (RyRy) Pn {acm -R.2)f 


+ 3Reo sin? oe : re \) 


Oo 


3. Ky = LL, - (Xn, + Xyo) 
Sharp Ogive (Figure I-2) 


1. Using parameters determined for volume of 


sharp ogive, calculate: 


= LL, - (2R° - 3 R,R* + R,°)/(3A,) 


sh bc Sen reel SS call 


deri SP sacra ae me als 


: 
i 
\ 
; 
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Spherically Blunted Ogive (Figure I-2) 


1. Using parameters determined for volume of 


blunt ogive, calculate: 


3 
LL 3 


rs 1 . 2. 3/2 _ 


Xn. = (same as for blunt cone) 
2. Xy = LL, - (Xy, + Xwo) 


Centroid of Body 


Cylinder 
.} 
rr B 
ne 2 
Centroid of Boattail 


Truncated Cone 
Xp = 


Centroid of Complete Body 


Lan fo ar+a 
[2 $9-+- 


a [ KpAag + (LL, + sts +(LL, + &, + Tan) Aga} 


P 


: enn see | So eT en 


~ SI EE ST ET TTI TT eR es eae eae oe a 3 ser 
Poa eT NH a — 


— 


aN A 


ve oe ET TE 
ee Na ea et aarti 


Blunt cone dimensions. 
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Figure I-2 e 
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NOMENCLATURE 


Fin forward wedge length, in. 
Plan area of body, ine” 
Plan area of boattail, aes 

Pla’ area of nose, in.? 

Plan area of sharp ogive, in.? 
Exposed aspect ratio = b2/ (2 S¢) 

Area involved in plan area calculation 
Area involved in plan area calculation 
Fin aft wedge length, in. 

Exposed fin semispan, in. 

Total span of fins plus body, in. 
Crossflow drag coefficient 

Lift coefficient 

Rolling moment coefficient 

Pitching moment coefficient 

Body alone pitching moment coefficient 
Hinge-moment coefficient 

Root bending-moment coefficient 

Normal force coefficient 

Body alone normal force coefficient 
Installed fin normal force coefficient 


Fin alone normal force coefficient 
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CNpp Fin in presence of body normal force coefficient 


Ca Yawing momern. ccefficient 

Cp Root crord, in. 

Cy. Vip chord, in. 

CP xuTA Fin alone centcr of pressure measured in X 
direction from hinge line, 1/C, 

CPXutp Fin in presence of body center of pressure 


measured 1n X direction from hinge line, 1/C, 
CPypca Fin alone center of pressure measured in Y 
direction from root chord, 1/(b/2) 
CPY RoR Fin in preserce of body center of pressure 


measured in Y¥ direction from root chord, 


1/ (b/2) 
a Body diameter, reference length, in. 
aT Base diameter, in. 
£ Body norma. force, lbs 
HL itinge line position from leading edge of root 


chord, 1/C,, 


Ry Apparent mass factor 

Ky Apparent mass tactor 

L. Distance from base of nose to origin of Ry in. 

L3 Distance from crigin ot Ry to theoretical nose, in. 
LI.) Distance from base of nose to perpendicular from 


tangent cf nose blunting te centerline, in. 


2 Actual nose length, in. 
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Total configuration length, in. 

Body length, in. 

Nose actual length, in. 

Nose theoretical length, in. 

Mach number 

Body pitching moment, in.-lbs 

Crossflow Mach number 

Calculated factor associated with nose planform 
area 

Calculated factor associated with nose planform 
area 

Tctal pressure, psfa 

Ogive radius, in. 

Crossflow Reynolds number 

Reynolds number based on bcay diameter 

Unit Reynolds number 

Ratio of R to LL, 

Ratio of Lup to LL, 

Ratio of L, to LL, 

Nose blunting radius, in. 

Ratio of Ry to LL, 

Ratio of R to LL, 


Distance from origin of R to origin of Rye in. 


Distance from origin of R to rose centerline, in. 
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XCP 


Cross sectional area of body, reference area, 
tae 

Area of base, in.” 

Exp2sed area of one fin, in.” 

Totai planform area, Swe? 

Thickness of fin at root, in. 

Thickness of fin at tip, in. 

Free stream velocity 

Total volume, in.? 

Volume of body, in.? 

Volume of boattail, in.? 

Volume of nose, in? 

Centroid of total plan area, measured from actual 
nose, in. 

Centroid of plan area of body, in. 

Centroid of plan area of boattail, in. 

Centroid cf plan area of nose, in. 

Centrcid of truncated nose, in. 

Centroid of nose blunting, in. 

Distance from nose to centr of pressure, calibers 

Distance fzom nose to body alone center of 
pressure, calibers 

Distance from moment reference to wffective center 


nf pressure of body on fin interference, calibers 
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XCPRrH 


Xcpra 


XCPrin 


XCPron 


Distance from hinge line to effective center of 
pressure of body on fin interfexence, i/c, 

Distance from moment reference to fin alone 
center of pvressure, calibers 

Distance from leading edge of root chord to 
center of pressure of installed fin, =/CR 

Distance from moment reference to etfective 
center of pressure of fin on body interference, 
calibers 

bistance from moment reference *o hinge line, 
negative aft of moment reference, calibers 

Distance from actual nose to moment reference, in. 

Distance from root chord to effective center of 
pressure of body on fin interference, 1/(b/2) 

Distance from root chord to center of pressure 
of installed fin, 1/(b/2) 

Ratio of Rp to LL, 

Empirical correction to body alcne pitching 
moment coefficient 

Maximum value of empirical picching moment 
correction 

Angle of attack 

Prebend angle 

Angle associated with nose plan area calculation 


Regression coefficient (intercept) 
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Regression voefficient (A) 

R2egresuion coefficient ay 

Regression coefticient (AR) 

Regression coefficient (d/b') 

Incremental normal force @ve to pedy on fin 
interference 

Snerymental april force due to fin on body 
interference 

Normalized hody alone pitching monent correction 

Ratts of Sroserice Grag of a circular cylinder of 
finite Jength to one Of infinite Jength 

Angle associatea with calculation of nose plan 
area, vadians 

Leading sdge sweep angle of fir, deg 

Taper ratio cr fin (Sp/Cp} 

Air density 

Standard deviation 

Roll angle, deg 

Angle associated with nose plan area calculation, 


radians 
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